JAGIELLONIAN UNIVERSITY 
INSTITUTE OF PHYSICS 




Search for the CV symmetry violation in the 
decays of Ks mesons using the KLOE detector 



Michat Silarski 



PhD thesis prepared in the Department of Nuclear Physics 
of the Jagiellonian University under supervision 
of Prof. Pawel Moskal 



Cracow 2012 



UNIWERSYTET JAGIELLONSKI 
INSTYTUT FIZYKI 




Poszukiwania lamania symetrii CP w rozpadach 
mezonu Ks za pomoc^ detektora KLOE 



Michal Silarski 



Praca doktorska wykonana w Zakladzie Fizyki J^drowej 
Uniwersytetu Jagiellonskiego pod kierunkiem 
Prof, dr hab. Pawla Moskala 



Krakow 2012 



„ Nur ein Lebenfur die anderen ist lebenswert " 

Albert Einstein 



Abstract 



The aim of this work was to determine the Ks Svr^ decay branching ratio and a 
modulus of the r/ooo parameter, defined as the ratio of amplitudes for Ks — )■ Sn^ to 
Kl — 37r° decays, which characterizes the CV symmetry violation in this decay. 
The measurement has been carried out with the KLOE detector operating at the - 
factory DA$NE in the Italian National Center for Nuclear Physics in Frascati. DA<I>NE 
collides the e"*" and e~ beams at the center of mass energy ^/s = 1019.45 MeV. The 
e+e~ collisions result in the meson creation which is almost at rest and decay pre- 
dominantly to kaon pairs. The decay products are registered using the KLOE detection 
setup, which consists of large cylindrical drift chamber surrounded by the electromag- 
netic calorimeter. The detectors are placed in a magnetic field of B ^ 0.52 T generated 
by superconducting solenoid. The Ks mesons were identified with high efficiency via 
registration of these Kl mesons which crossed the drift chamber without decaying and 
then interacted with the KLOE electromagnetic calorimeter. The Ks four - momen- 
tum vector was then determined using the registered position of the K^ meson and the 
known momentum of the meson. Next, the search for the Ks 37r° — t- 67 decay 
was carried out by the selection of events with six gamma quanta which momenta were 
reconstructed using time and energy measured by the electromagnetic calorimeter. To 
increase the signal over background ratio after identification of the Ks meson and re- 
quiring six reconstructed photons a discriminant analysis is performed. It is based on 
kinematical fit, testing of the signal and background hypotheses and exploiting of the 
differences in kinematics of the Ks decays into 27r° and 37r°. 

The search for the Ks 37r° decay presented in this work failed to detect a signal 
of sufficient statistical significance. Hence, we have obtained the upper limit on the 
Ks 37r° branching ratio at the 90% confidence level: 

BR{Ks 37r°) < 2.7 ■ 10^*^ , 

which is almost five times lower than the latest published result. This upper limit 
can be translated into a limit on the modulus of the 77000 parameter amounting to: 
|'7ooo| < 0.009 at the 90% confidence level. This corresponds to an improvement of 
the |?7ooo| uncertainty by a factor of two with respect to the latest direct measurement. 
The upper limit on the Ks 37r° branching ratio determined in this work is still 
about one order of magnitude larger than the prediction based on the Standard Model. 
Hence, the search for this decay will be continued with the upgraded KLOE detector, 
which has realistic chances to observe the Ks 37r° decay for the first time in the 
near future. 



Streszczenie 



Celem tej pracy bylo wyznaczenie stosunku rozgalgzieri dla rozpadu Ks — )■ 37r°, oraz 
modulu parametru tjqqq, zdefiniowanego jako stosunek amplitud na rozpad Ks — )■ 37r° 
i Kl — 37r°, ktory charakteryzuje niezachowanie symetrii CV w tym procesie. 
Pomiary wspomnianego rozpadu byly prowadzone za pomocq detektora KLOE dziala- 
j^cego na akceleratorze wi^zek przeciwbieznych DA$NE we Wloskim Narodowym 
Centrum Fizyki Jqdrowej we Frascati. DA$NE zderza wiqzki e" przy energii 
w centrum masy ^/s = 1019.45 MeV rownej masie mezonu 0. W wyniku zderzen 
e~^e~ powstaj^ mezony (p. Produkowane one praktycznie w spoczynku i rozpadaj^ 
si? glownie na pary kaonow. Do pomiaru powstaj^cych czqstek wykorzystywany jest 
uklad detekcyjny KLOE. Zbudowany jest on z cylindrycznej komory dryfowej otoczo- 
nej kalorymetrem elektromagnetycznym. Calosc umieszczona jest w polu magnety- 
cznym nadprzewodzqcego solenoidu o indukcji B 0.52 T. Mezony Ks identyfi- 
kowane byly z duz^ wydajnosciq poprzez rejestracj? tych mezonow Kl, ktore nie 
rozpadly sig w komorze dryfowej i zareagowaly z materialem kalorymetru elektro- 
magnetycznego. Wektor czteropgdu mezonu Ks okreslany byl na podstawie zareje- 
strowanej w kalorymetrze pozycji mezonu Kl i znanego wektora czteropgdu mezonu 
(p. Poszukiwania rozpadow Ks — )■ 3vr° — )■ 67 prowadzone byly nastgpnie poprzez wy- 
branie zdarzen z szescioma zrekonstruowanymi kwantami gamma. Ich p?dy okreslane 
byly na podstawie czasu i energii mierzonych za pomocq kalorymetru. Aby zwigk- 
szyc stosunek sygnalu do tla zdarzenia ze zidentyfikowanym mezonem Ks i szes- 
cioma zrekonstruowanymi kwantami gamma zostaly poddane dalszej analizie, opartej 
na dopasowaniu kinematycznym, testowaniu hipotez tla i sygnalu oraz wykorzystaniu 
roznic w kinematyce rozpadow mezonu Ks na dwa i trzy mezony 7r°. 
W wyniku przeprowadzonych poszukiwari nie zarejestrowano zadnego zdarzenia od- 
powiadajqcego rozpadowi Ks 37r°. Dlatego okreslono gom^ granicg stosunku roz- 
galgzien dla rozpadu Ks — )■ 37r° na poziomie ufnosci 90%: 

BR{Ks 37r°) < 2.7-10^*^ . 

Otrzymana wartosc jest okolo pigc razy nizsza od ostatniej opublikowanej gomej gra- 
nicy na ten stosunek rozgalgzieri. Otrzymana na poziomie ufnosci 90% wartosc gomej 
granicy dla modulu parametru r/ooo wynosi: |?7ooo| < 0.009, co odpowiada zmniejsze- 
niu jego niepewnosci dwa razy w stosunku do ostatniego bezposredniego pomiaru. 
Otrzymana wartosc gomej granicy na BR(Ks 37r°) jest wigksza od teoretycznych 
przewidy wan o rz^d wielkosci, dlatego poszukiwania tego rozpadu bgdq kontynuowane 
za pomocq systemu detekcyjnego KLOE wyposazonego w nowe detektory. 
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Chapter 1 



Introduction 



In 1918 Amalie Emmy Noether proved a theorem connecting the symmetries of the 
physical systems and conservation laws [1] which became one of the greatest achieve- 
ments of the twentieth century theoretical physics. It shows for instance that a system 
invariant under translations of time, space, or rotation will obey the laws of conserva- 
tion of energy, linear momentum, or angular momentum, respectively. From that time 
on the symmetries have become an essential part of almost all physics theories and 
models, especially in the particle physics. And so for example every quantum field 
theory describing the interaction and properties of elementary particles are formulated 
requiring the Lorentz invariance. Furthermore, the discrete symmetries of Parity V, 
Charge Conjugation C and Time Reversal T proved to be very useful in the calcula- 
tion of the cross sections and decay rates, especially for the processes governed by the 
strong interaction. These symmetries became also an important tool in the Standard 
Model formulation. 

Among the known elementary forces the weak interaction has appeared to be very 
peculiar mainly because it violates V and C symmetries [2, 3] as well as their combi- 
nation: CV. The CV violation was discovered unexpectedly in 1964 by Christenson, 
Cronin, Fitch and Turlay during the regeneration studies of the neutral K mesons [4]. 
In the framework of Standard Model the CV violation mechanism is introduced by the 
quark mixing described by the complex Cabibbo - Kobayashi - Maskawa matrix with 
one nonzero phase [5, 6], which requires the existence of three generation of quarks. 
Parameters describing the CV violation in the neutral kaon system were measured with 
a good precision by several experiments [7-9], and at present the main experimental 
effort is focused on studies of the neutral B and D meson systems [10-13]. However, 
there are still several interesting open issues in the kaon physics. One of them is the 
Kg Stt" decay which, assuming the CVT invariance, allows one to investigate the 
direct CV symmetry violation [14]. Despite several direct searches [15, 16] and KsKl 
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interference studies [17, 18], this decay remains undiscovered and the best upper limit 
on the branching ratio BR{Ks Svr") < 1.2 ■ 10~'^ [16, 19] is still two orders of 
magnitude larger than the predictions based on the Standard Model [14]. 

This work is focused on the measurement of the Ks — )■ Svr*^ decay branching ratio 
based on the data sample gathered in 2004 - 2005 with the KLOE detector operat- 
ing at the (j) - factory DA$NE in the Italian National Center for Nuclear Physics in 
Frascati. DA$NE collides the e"*" and e~ beams at the center of mass energy of -y/i 
= 1019.45 MeV near the (j) meson mass [20]. The e+e~ collisions result in cf) meson 
creation which is almost at rest {(3^ ^ 0.015) and decay predominantly to K^K~ 
(49%), KsKl (34%), pvr (15%) and r/7 (1.3%) final states [19]. The decay products 
are registered using the KLOE detection setup, which consists of large cylindrical drift 
chamber surrounded by the electromagnetic calorimeter. The detectors are placed in 
a magnetic field of i? ~ 0.52 T generated by superconducting solenoid. Since the (p 
mesons are produced almost at rest, kaons arising from the decay move with the rel- 
ative angle close to 180°, and as a consequence, their decay products are registered 
in the well separated parts of the detector. The Kg mesons are identified with high 
efficiency (~ 34%) via registration of these A'l mesons which cross the drift cham- 
ber without decaying and then interact with the KLOE electromagnetic calorimeter (so 
called Ks tag). The Kg four - momentum vector is then determined using the reg- 
istered position of the K^ meson and the known momentum of the meson, which 
is estimated as an average of the momentum distribution measured using large angle 
e+e" scattering. The search for the Kg — )■ 37r° — )• 67 decay is then carried out by 
the selection of events with six 7 quanta which momenta are reconstructed using time 
and energy measured by the electromagnetic calorimeter. Background for the searched 
decay originates mainly from the Kg — )■ 2n^ events with two spurious clusters from 
fragmentation of the electromagnetic showers (so called splitting) or accidental activ- 
ity, or from false Kl identification [21]. To increase the signal over background ratio 
after identification of the Kg meson and requiring six reconstructed photons a dis- 
criminant analysis is performed. It is based on kinematical fit, testing of the signal and 
background hypotheses and exploiting of the differences in kinematics of the Kg de- 
cays into 27r° and 37r°. 

This thesis is divided into nine chapters. The detailed description of the CV viola- 
tion mechanism in the neutral kaon system is presented in the second chapter together 
with the motivation to search for the Kg — > 37r° decay. 

The third chapter provides the description of experimental tools used for the measure- 
ment: the DA<I>NE collider, the KLOE detector as well as the trigger and data acquisi- 
tion system. 
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The method used to identify the Ks mesons based on the detection of the interac- 
tions in the electromagnetic calorimeter is presented in chapter four. 
Chapter five is devoted to the identification of the Ks 27r° events used for the nor- 
malization of the measured branching ratio. 

The discriminant analysis used to increase the signal over background ratio is described 
in chapter six, where also the background estimation based on the Monte Carlo simu- 
lations and the final result of the Ks — )■ 37r° identification are presented. 
The seventh chapter is devoted to the estimation of the systematic uncertainties of the 
measurement. 

The determined upper limit on the Ks — )■ 37r° branching ratio is given in chapter eight 
together with the estimation of the modulus of the r^ooo parameter, defined as the ratio 
of amplitudes for Ks — Stt" to Kl 37r° decays. 

Finally, the ninth chapter comprises the summary and perspectives. In particular we 
discuss the possibility of the first observation of the Ks Stt^ decay in the next 
KLOE-2 data - taking campaign during 2013 - 15. It will be conducted with the KLOE 
detector upgraded by Inner Tracker and with improved photon acceptance brought 
about by new calorimeters installed in the final focusing region. 



Introduction 



Chapter 2 

Neutral Kaon system and CV violation 



Discrete symmetries as parity V, charge conjugation C and time reversal T, as well as 
their combinations such as CV and CVT play a fundamental role in particle physics. 
The parity transformation V changes the signs of the three space coordinates, while C 
changes particle to its antiparticle and vice versa changing its internal quantum num- 
bers. The strong and electromagnetic interaction preserve eigenvalues of both, parity 
and charge conjugation operators, as well as eigenvalues of the CV operator. The weak 
interaction instead do not preserve these quantum numbers which result in a far - 
reaching consequences, especially in case of the CV operator. In the framework of 
the Standard Model the CV violation implies the existence of the third generation of 
quarks. Moreover, it is a very important mechanism which could have an essential con- 
tribution to the asymmetry between matter and antimatter in the Universe. In 1967 A. 
Sakharov laid out three conditions that would enable a Universe containing initially 
equal amounts of matter and antimatter to evolve into a matter dominated universe, 
which we see today [22]. The first condition was a violation of the baryon number 
conservation B, for which there is still no experimental evidence. However, simple 
baryon number violation would not be enough to explain matter - antimatter asymme- 
try if C and CV were exact symmetries ^ As the third condition Sakharov proposed that 
the Universe should undergo a phase of extremely rapid expansion [22]. 
Since the first discovery of the CV - violating neutral kaon decay in 1964, there have 
been made a big effort to describe the CV symmetry breaking within the Standard 
Model. The favoured theoretical framework was provided in 1 973 by Kobayashi and 
Maskawa, who pointed out that CV violation would follow automatically if there were 
at least six quark flavours. At present the main experimental effort is focused on the 
neutral B and D meson system studies [10-13]. However, there are still several inter- 
esting open issuses in the kaon physics which, as it will be shown in this chapter, can 



'in that case there would be a perfect equality between rates of different B violating processes and 
no asymmetry could be generated from the initially symmetric state [23]. 
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contribute to our better understanding of the CV violation mechanism. 



2.1 The neutral kaon system 

Kaons were discovered in 1947 by G. D. Rochester and Clifford C. Butler while study- 
ing cosmic ray showers with a cloud chamber [24]. The contrast between the produc- 
tion and decay times of these new particles entailed in 1953 introduction of a new 
quantum number called „strangeness" S [25, 26]. Strangeness is conserved by both 
electromagnetic and strong processes while first order weak interaction can induce 
transitions with AS = 1. 

Kaons appear in isospin I = \ doublets: (A'+, A'°) with S = \ and (l^ ,K~) with 
S = —1. They can be produced via strong interactions in processes e.g. like: 

vr+p — )■ K^K^p 

pp K"TT+K° 
PP — )■ f^+TT" A^ . 

From the point of view of strong interactions the meson is a particle with a cor- 
responding antiparticle K*^ . Violation of strangeness conservation by weak interaction 
allows for transitions like — t- 27r — )• A^ or K'^ — Stt — Ji°. Thus, the two 
strangeness eigenstates can oscillate one into another via the AS = 2, second order 
weak interactions, i.e., via virtual 2-pion and 3-pion states. The corresponding quark 
diagrams for these transitions are presented in Fig. 2.1. Therefore, in the evolution of 
a kaon in a free space states with well defined mass and width are mixtures of and 

[27]. The time evolution of the neutral kaon system, which may be given in the 
- rest frame is determined by Hamiltonian H and the following equation: 



M - ^- , (2.1) 






where M and F are 2 x 2 hermitian mass and decay matrices, respectively. In the Weis- 
skopf - Wigner approximation the elements of the mass matrix M can be expressed as 
a sum of contributions due to strong and weak interactions [27]: 

Mij = rrikdij + Hw j) + > . F ' (^.2) 
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Figure 2.1: Quark diagrams for the K - K transitions 



where is the neutral kaon mass and Hw denotes the effective weak Hamiltonian. 
The sum J2n 1^) (^1 ^^^^ virtual and real states connecting and A'°. The de- 
cay matrix elements related to the kaon decay width by unitarity, originate instead 
only from Hw [28]: 



Tij = 2n ^ 5{En - mA')(^|Hw|?^)(?^|Hw|j) • 

n^KO.K° 



(2.3) 



If the Hamiltonian of the system is invariant under T, CV and CVT Mij and Ty have 
to satisfy the following relations: 

|Mi2 - 2Fi2/2| = \M^2 - ^r];2/2| (T conservation), 

|Mi2 - iFi2/2| = |Mi*2 - i^\2l'A and Mn = M22, Fn = F22 {CV invariance) 
Mil = ^22; Til = F22 (CPT conservation). 

Without any assumption about symmetry invariance the Hamiltonian eigenstates of the 
neutral kaon system seen in nature can be written in the following form [29] : 



1 



V'2(rTN? 
1 

www 



;i + 6^)1/0 + (1- 65) liT) 



;i + e,)|ir°)-(i-6,)|r) 



(2.4) 



where eg and ei, are complex parameters expressing possible CV and CVT violation. 
In particular, if CVT invariance holds: = ez, = e. It is important to stress, that Ks 
and A'i are kaon states which preserve their identity during the evolution in free space. 
The experimental values for the lifetimes of these two particles differ by three orders 
of magnitude. The lifetime of the „short" state IA'5) amounts to 
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Ts = (8.953 ± 0.005) ■ 10"^^ s, while the „long" living particle \Kl) has a lifetime 
tl = (5.116 ± 0.020) ■ 10~^ s [19]. This large difference was explained by assuming 
CV to be an exact symmetry of the weak interactions. In this case the mass eigenstates 
defined in Eq. 2.4 reduce to the CV eigenstates^ (e^ = = 0): 



\Ki) 



V2 



\K') + W) 

|A'0) - 1:^) 



with CV = 1 
with CV = -1. 



(2.5) 



Neutral kaons decay mainly to the two - and three - pion final states with a well defined 
CP eigenvalues [30]: 

CVln+n-) = \n+n-) 

0^0\ _ |^0^0\ 



CV 



n IT TT 



') = (-1) 



l+l 



TT ' TT 



CP|7r"7r"7r") = -|7r"7r"7r") . 

For the Itt+tt^tt^) final state the eigenvalue depends on the total angular momentum 
/. However, since the three pions from the kaon decay are mainly in the relative s - 
wave state we can assume with a good approximation that the (tt'^, tt", 7r°) system is 
CV - odd. Thus, CV conservation would imply that | A'l) state is allowed to decay only 
to two pions while the „long" living IA'2) decays only to three pions state. Moreover, 
the large phase space difference between these two decay modes manifests itself in the 
difference between observed lifetimes. 

However, as it is presented in the next section, the CV invariance is violated by the 
weak interaction which entails big consequences for the whole particle physics and 
cosmology. 



2.2 CV violation in kaon decays 

In 1964 an experiment by Christenson, Cronin, Fitch and Turlay, unexpectedly ex- 
hibited that the long - lived kaon can decay also to the two - pion final states with 
branching ratio of about 2 ■ 10^^ [4]. Thus, the neutral kaons states seen in nature are 
not CP eigenstates defined in Eq. 2.5. However, they still can be expressed in the (| A'l), 

2 Here we assume a phase convention where CVlK^) = \K") and CVlK^) = \K°) [29]. 
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\K2)) basis: 



1+ e 



2 



\Ks) = -7^^=^ {\Ki) - e\K2)) . (2.6) 



1+ e 



Since up to now there is no signs of the CVT symmetry violation from now on we 
assume ^ eg = = e. 

We can understand the CV symmetry breaking within the scope of two distinct mech- 
anisms referred to as „direct" and „indirect" breaking. The „indirect" violation corre- 
sponds to the statement that the true eigenstates of both the strong and electroweak 
interactions are not exactly CV eigenstates but have small admixtures of the state with 
opposite CV [30]. It is also possible that CV violation occurs „directly" in the weak 
decays themselves. In the following the „direct" violation will be explained on the ex- 
ample of kaon decays to two pions. The two - pion systems originating from decays 
of kaons can be produced with isospin / = 0or/ = 2( isospin equal one is forbidden 
by Bose symmetry [30]): 

k°vr°) = 1 = 0)- ^Jlw'A I = 2) 

Itt+tt-) = y||vr+7r-; / = 0) + y||vr+7r-; 1 = 2). 

The corresponding weak decay amplitudes of kaons can be expressed as [28]: 

(7r7r;/|Hw|i^°) = Aje'^' 
(7r7r;/|Hw|^) = A*je'^' , 

where we have explicitly exhibited the final state phases 5[, which arise from the fi- 
nal state strong interactions of the pions. Direct CV violation, occurring at the decay 
vertices, appears as a phase difference between the Aq and A2 amplitudes'^ [28]. This 
phase difference is generated by a class of so called „pengiun" diagrams for s quark 
decay, one of which is presented in Fig. 2.2. 

Typically the CV violation in the neutral kaon sector is characterized in terms of the 
following parameters: 

A(Kl TT + TT") , , ^ 

= a)t^^ = ^00 e^- = e - 2e', (2.7) 



^Although there are some theoretical predictions for the CVT violation [31-34], all the tests done 
so far resulted in the confirmation that it is an exact symmetry [19, 35-38]. 
"^In general ^0 and ^2 could be real if CV would be conserved. 
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Figure 2.2: „Pengiun" diagram for the s quark decay. 



where e is the mixing parameter defined before in Eq. 2.6 and e accounts for the direct 
CV violation and can be expressed in terms of the weak amplitudes [28]: 

, (7r7r;0|Hw|/^5)(7rvr;2|Hw|A^L) - (tttt; 0|Hw|i^L) (vrvr; 2|Hw|i^5) 



v^(7r7r;0|Hw|/^5)' 



?a i ■= — Im — — 

y/2 \Aq 



(2.8) 



Above defined parameters were measured many times and are known with a good 
precision [19]: 



|r/+_| = (2.232 ±0.011) ■ 10 
|r/oo| = (2.221 ±0.011) ■ 10 
|e| = (2.228 ±0.011) • 10 



^_ = (43.51 ±0.05)° 
(,00 = (43.52 ± 0.05)° (2.9) 
= (43.51 ±0.05)°. 



Moreover, measurements of the double ratio of the two pion decay rates |77+_|/|77oo| 
have proved that e' is different from zero indicating occurrence of the direct CV viola- 
tion [38]. 

Analogous CV invariance breaking should appear also in the Ks decays. As before we 
can define the following amplitude ratios; 

A(Ks vt+tt-ttO) , , ^ 



_ A{Ks ^ nW) _ ,,000 ^,,,/ .2 11) 



As in the case of two - pion decays the ratios contain direct CV violation parameters 
related in the lowest order of the Chiral Perturbation Theory by the following equa- 
tions: e'^_o = ^000 ~ ~2e' [39]. The possible |7r+7r~7r°) final state originating from 
the neutral kaon decays can be produced with isospin / = 0, 1, 2, or 3. The / = 
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and 1 = 2 states have CV = 1, and Ks can decay into them without violation of the 
CV symmetry. However, they are expected to be strongly suppressed by centrifugal 
barrier effects [19]. For the 1=1 and / = 3 states there is no centrifugal barrier and 
CV = — 1 so Ks decay requires violation of this symmetry. Anyhow the two kinds of 
final states can be separated by the analysis of the 7r+7r~7r° Dalitz plot [19]. In the case 
of |7r°7r°7r°) final state, only isospin / = 1 or / = 3 is allowed, for which CV = — 1. 
Therefore, the Ks — )■ 37r° decay is a purely CV violating process [19]. 
The present knowledge about ?7+_o and ?7ooo is poor mainly due to very low decay rates 
for the Ks 37r decays. The current value of the Ks — ?■ 7r+7r~7r° branching ratio 
amounts to BR{Ks n+n-n^) = {3.5tol) ■ 10"^, and the Ks 37r° has been 
never observed [19]. The best upper limit on this decay branching ratio was set by 
KLOE collaboration and amounts to BR{Ks 37r°) < 1.2 ■ 10"^ [16], while the 
prediction based on Standard Model is equal to about 2 ■ 10"^ [39]. The corresponding 
knowledge about the amplitude ratios can be summarized as follows [19]: 

i?e(r/+_o) = -0.002 ± 0.007+°;°°^ 
/m(r/+_o) = -0.002 ± 0.009 
Im{r]ooo) = -0.001 ± 0.016 
|r/ooo| <0.018. 

Therefore, it is clear, that the full understanding of the CV violation in the neutral kaon 
system demands new high statistics measurements, in particular of the Ks — )■ 37r° 
decay which is a subject of this work. One of the high precision experiments, which 
has been greatly contributed to this quest, is the KLOE detection setup which will be 
presented in the next chapter. 
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Chapter 3 

The KLOE experiment at DA$NE 



In this chapter the characteristics of the DA$NE collider and the KLOE detector are 
briefly described. More detailed description can be found in Ref. [20, 40-44]. 

3.1 The DA<I>NE Collider 

DA$NE is an e'^e~ collider, optimized to work with a center of mass energy around 
the (j) mass, = (1019.418 ± 0.008) MeV [43]. The „heart" of the collider are two 
storage rings in which 120 bunches of both, electrons and positrons, are stored. Each 
bunch collides with its counterpart once per turn, minimizing the mutual perturbations 
of colliding beams. Electrons are accelerated to final energy in the Linac (see Fig. 3. 1), 




Figure 3.1: Schematic view of the DA<I>NE collider. The picture was adapted from [20]. 

accumulated and cooled in the accumulator and transferred to a single bunch in the 
ring [20]. Positrons are created in an intermediate station in the Linac, and then follow 
the same procedure as electrons. Both, electrons and positrons are injected into the 
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rings at final energy of about 510 MeV. The beams collide in two interactions regions, 
with a frequency up to 356 MHz, corresponding to a minimum bunch crossing period 
of Trf = 2.7 ns. The e^e" collisions result in meson creation which is almost at 
rest {(3^ ^ 0.015) and decay predominantly to K+R- (49%), KsKl (34%), pn (15%) 
and ?77 (1.3%) final states [19]. The decay products are registered using the KLOE 
detection setup, which consists of large cylindrical drift chamber surrounded by the 
electromagnetic calorimeter. The components of KLOE will be briefly described in 
the next section. 



3.2 The KLOE detector 

The KLOE detector consists of a large cylindrical drift chamber and a hermetic electro- 
magnetic calorimeter. A superconducting coil and an iron yoke (see Fig. 3.2) surround- 
ing the calorimeter provides a 0.52 T magnetic field. The beam pipe at the interaction 
region is a beryllium sphere with 10 cm of radius and 0.5 mm thick. This structure min- 
imizes both the multiple scattering and the energy loss of the charged particles from 
Ks decays, as well as the probability of Kl regeneration [43]. 





Figure 3.2: Vertical cross section of the KLOE detector. The figure is adapted 
from [20]. 
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3.2.1 The Drift Chamber 

The KLOE drift chamber has a cylindrical shape 3.3 m long, with internal and exter- 
nal radii of 25 cm and 2 m, respectively [43]. It was designed to register all charged 
secondary products from the Kl decay and measure their properties with great preci- 
sion [20]. Thus, its size was dictated by a long lifetime of this particle^ To minimize 
the Kl regeneration, multiple Coulomb scattering and photon absorption KLOE drift 
chamber is constructed out of carbon fiber composite with low-Z and low density, and 
uses a gas mixture of helium (90%) and isobutane (10%) [21]. The radiation length of 
the gas amounts to about 900 m, including the contribution of the 52140 wires [43]. 
In order to obtain high and uniform track and vertex reconstruction efficiencies, wires 
are strung in an all - stereo geometry, with stereo angles varying with the radius from 
50 mrad to 120 mrad going outward [43]. This design results in a uniform filling of 
the sensitive volume with almost square drift cells, with shape slowly changing along 
z axis^. Fig. 3.3 shows the wire geometry during the drift chamber construction as il- 
luminated by light. Particles from the decays are produced with small momenta and 
therefore track density is much higher at small radii [43]. Thus, dimensions of the cells 
were designed to be of about 2x2 cm^ for the 12 innermost wire layers, and to of 
about 3x3 cm^ for the remaining 48 layers [42]. 

To extract the space position from the measured drift time of the incident particle, 232 
space - to - time relations are used. They are parametrized in terms of two angles /3 
and (j) defined in Fig. 3.4. The /3 angle characterizes the geometry of the cell directly 
related to the electric field responsible for the avalanche multiplication mechanism, (j) 
instead gives the orientation of the particle trajectory in the cell's reference frame, de- 
fined in the transverse plane and with origin in the sense wire of the cell [20]. 
Using the wire geometry, space - to - time relations and known magnetic field one 
can reconstruct the tracks and vertices of charged particles. The reconstruction pro- 
cedure starts with pattern recognition and is followed by track and vertex fitting. The 
pattern recognition associates hits^ close in space to form track candidates and gives 
a first estimate of the track parameters. Then track fitting provides the final values of 
these parameters by minimization procedure based on the difference between the fitted 
and the expected drift distances (so called residuals), as evaluated from measured drift 
times and space - to - time relations. Finally the vertex fit procedures search for pos- 
sible primary and secondary vertices, on the basis of the distance of closest approach 
between tracks [20]. 

To ensure the stability in time of the KLOE drift chamber performance, the system 

'The mean decay path of the meson produced in the decay amounts to about 3.4 m [20]. 
^ The z axis of the KLOE reference frame is defined as the bisector of the angle between colHding 
e+ and beams [20]. 

^As a hit we consider a presence of a signal on a sense wire. 
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Figure 3.3: Drift chamber stereo wires geometry. The figure is adapted from [45]. 




Figure 3.4: Left: Wire geometry with the definition of stereo angle e between the wire 
of length L and the z - axis. Right: Definition of (3 and angles characterizing the 
shape of the cell and the angle of the incident track [20]. The figures are adapted 
from [20]. 
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is calibrated periodically by acquiring samples of cosmic ray events suitable for the 
measurement of about 200 different space - to - time relations [20]. The calibration 
is performed at the beginning of each KLOE run and selects about 80000 cosmic ray 
events [42] . These events are tracked using the existing space - to - time relations and 
the average value of the residuals for hits in the central part of the cells is monitored. If 
the residuals exceed 40 i^im additional 3-10^ cosmic ray events are collected, and a new 
set of calibration constants is obtained. Finally, during data taking the drift chamber 
performances are monitored using selected samples of events [42]. 
The KLOE drift chamber provides tracking in three dimensions with a resolution in 
the transverse plane of about 200 fim, resolution in the z-coordinate measurement of 
about 2 mm and of 1 mm on the decay vertex position. The momentum of the particle 
is determined from the curvature of its trajectory in the magnetic field with a fractional 
accuracy ap/p = 0.4% for polar angles larger than 45° [20]. 

3.2.2 The Electromagnetic Calorimeter 

The KLOE electromagnetic calorimeter was designed to provide hermetic detection of 
low energy gamma quanta with high efficiency, good energy resolution and excellent 
time resolution for the neutral vertex reconstruction and to trigger the events [20]. It 
consists of a barrel built out of 24 trapezoidal shaped modules and side detectors (so 
called endcaps) read out from both sides by a set of photomultipliers. The barrel is a 
cylinder with an inner diameter of 4 m, made of 24 modules 4.3 m long and 23 cm 
thick. Each endcap consists of 32 vertical C - shaped modules. This structure covers 
98% of the full solid angle. Each module consists of a mixture of lead (48% of the 
volume), scintillating fibers (42%), and glue (10%) [43]. Fibers, each with a diameter 
of 1 mm, are embedded in 0.5 mm lead foils accelerating the showering processes. The 
special care in design and assembly of the Pb - fiber composite ensures that the light 
propagates along the fiber in a single mode with velocity ~17 cm/ns, which greatly 
reduces spread of the light arrival time at the fiber ends [20]. Calorimeter modules 
are read out at both ends viewed by light guides of area of 4.4 x 4.4 cm^ coupled 
to the photomultipliers transforming the light into electric impulses. This defines so 
called „calorimeter cells" which form five larger structures (see Fig. 3.5): planes 4.4 cm 
wide"*. 

When a particle hits the calorimeter for each cell both, the charge as well as time of 
arrival of the photomultiplier signals are registered. The cell energy is taken as the 
average of the energy registered at both sides, after correcting for the light attenuation 
along the fiber [43]. The energy calibration starts by a first equalization in cell response 
to minimum ionizing particles at calorimeter center, and by determining the attenuation 

"*The last plane of cells is 5.2 cm wide. 
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Figure 3.5: Schematic view of the readout structure on one side of the barrel mod- 
ule [46]. 60 defined cells form 5 planes and 12 columns of the calorimeter module. 
Filled circles represents photomultipliers. The figure is adapted from [46]. 




Figure 3.6: Photograph of the KLOE calorimeter. One can see 24 modules of the barrel 
and the inner plane of one of the endcaps. The figure is adapted from [46]. 

length of each single cell using cosmic rays acquired in dedicated runs. This is done 
before the start of each long data taking period [20]. The energy determined from the 
measured amplitudes of signals for both sides A and S of a cell S"^^ amounts to: 

CA,B _ cAB 

E^'^(MeV) = — ^ ■ K, (3.1) 

where Sq ' are the zero - offsets of the amplitudes scale, Sm corresponds to the re- 
sponse for the minimum ionizing particle crossing the calorimeter center and K factor 
gives the energy scale in MeV [46]. The total energy deposited in a cell is calculated 
as the mean of values determined at both ends for each cell. The determination of the 
absolute energy scale K relies on a monochromatic source of 510 MeV photons: the 
e"*"e^ —7- 77 sample. The latter calibration is routinely carried out each 200 - 400 nb~^ 
of collected luminosity [20]. 
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For each cell two time signals T"^ and (digitized by the Time to Digital Con- 
verter TDC ) are recorded. The arrival time t and position s of the impact point along 
the fiber direction can be determined as^: 

t(ns) = l(t^ + t^-t^-t^)-|-, (3.2) 
.(cm) = |(t^-t^-t^ + t^), (3.3) 

with t^'^ = c^'^ • T^'^, where c^'^ are the TDC calibration constants, t^'^ denotes 
overall time offsets, L stands for length of the cell (cm) and v is the light velocity in 
fibers (cm/ns) [40]. 

Based on the reconstructed energies, times and positions cells are merged into calorime- 
ter clusters. First the adjacent cells are grouped into so called „preclusters"^. The time 
spread of cells forming the precluster has to be smaller than 2.5 ns [46] . Moreover, cells 
are merged in one cluster if a distance between them and the center of the precluster is 
less than 20 cm [46]. The cluster energy is evaluated as the sum of the cells energies: 



E,i = Y,Ei, (3.4) 



while the time and position centroids are obtained as weighted averages: 

= ^if'' (3.5) 

Rci = ^^V^ ■ (3.6) 

i denotes the ith cell belonging to the cluster and rj stands for the cell's position vector 
with respect to the interaction point. is next related to the time of flight of parti- 
cle from the interaction point to the cluster position. It is done subtracting the event 
global time offset, common to all channels and depending on the trigger formation 
time with respect to the real e+e~ interaction time. Due to the spread of the particle's 
arrival times, the KLOE trigger is not able to identify the bunch crossing related to 
each event, which has to be determined offline [20]. The common „Start" signal to the 
calorimeter TDC boards is provided by the first level trigger, which will be described 
in the next section. The „Stop" instead is given by the photomultiplier signals delayed 
because of the electronics and light propagation in the fibers [43]. Time measured by 
the calorimeter can be expressed as: 

Td = Ttof + 5c- Nbc ■ Trf , (3.7) 



''A cell is added into a precluster only if its times and energies were reconstructed at both sides of 
the calorimeter module. 
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where Ttof is the time of flight of a particle from the interaction point to the cluster 
position, 5c is a single number accounting for the overall electronic offsets and cable 
delays, and Ni^c is the number of bunch - crossing periods needed to generate the 
TDC start. The values of 5c and Ni^c are determined for each data taking run with 
e+e~ — 77 events by looking at the ^tof = Td — Rd/c distributions (c denotes 
the speed of light ) [20]. For such events this distribution shows well separated peaks 
corresponding to different values of Ni,c- We define 5c as the position of the largest 
peak in the distribution, and obtain T^/ from the distance between peaks [44]. During 
offline processing, to allow the cluster times to be related to the particle time of flight, 
we determine for each event the corrected cluster times: 

tci = Tel - {5c - Nhc ■ Trf) . (3.8) 



The KLOE electromagnetic calorimeter allows for measurements of particle energies 
and flight time with accuracies of as = -^^^=E and ait) = /"^p" © 140 ps, 

W E [GcV] W E [Ge V] 

respectively [40]. Analysis of the signal amplitude distributions allows to determine 
the location where the particle hits the calorimeter module with accuracy of about 1 
cm in the plane transverse to the fiber direction. The longitudinal coordinate precision 
is energy dependent: = ^'^^^^^ [21]. 

3.2.3 The Trigger system 

The KLOE trigger system is based on local energy deposits in the electromagnetic 
calorimeter and hit multiplicity information from the drift chamber. It has been opti- 
mized to retain almost all e+e" — )■ (j) decays, and provide efficient rejection on the 
two main sources of background: small angle e+e~ — > e^e~ scattering and particle 
lost from the DA$NE beams [20]. Moreover, all e"'"e~ scattering and 77 events pro- 
duced at large polar angles are gathered for detector monitoring and calibration. Since 
the DA$NE bunch crossing period amounts to T^f = 2.7 ns, KLOE trigger must op- 
erate in continuous mode. A two level scheme was chosen. A first level trigger Tl 
is produced with a minimal delay (~ 200 ns) and is synchronized with the DA$NE 
master clock [41]. The Tl signal initiates conversion in the front - end electronics 
modules, which are subsequently read out following a fixed time interval of about 2.6 
/iS. This corresponds to the typical drift time of electrons travelling in the drift chamber 
cells [20]. After the arrival of a first level trigger, additional information is collected 
from the drift chamber, which is used together with the calorimeter information as a 
second level trigger T2. It confirms the first level trigger, initializes digitisation of the 
drift chamber electronics and starts the data acquisition readout. If no T2 signal arrives 
before the end of 2.6 yus dead time, all readout is reset [41]. 
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Tl and T2 triggers are based on the topology of energy deposits in the KLOE electro- 
magnetic calorimeter and on the number and spatial distribution of the drift chamber 
hits. Since decay events have a relatively high multiplicity, they can be efficiently 
selected by the calorimeter trigger by requiring two isolated energy deposits above a 
threshold of 50 MeV in the barrel and 150 MeV in the endcaps. Events with only two 
fired sectors in the same endcap are rejected, because this topology is dominated by 
machine background. Moreover, we require about 15 hits in the drift chamber within a 
time window of 250 ns from beam crossing [20]. The trigger identifies e^e~ — )■ e^e~ 
events requiring clusters with energy of about 350 MeV. An event which satisfies at 
least one of the two above conditions and is not recognized as e'^e~ scattering, gen- 
erates a first level trigger Tl^. The level - 2 trigger T2, requires further multiplicity 
or geometrical conditions for the electromagnetic energy deposits, or about 120 drift 
chamber wire signals within a 1.2 /is time window. At the level 2 trigger recognizes 
also the cosmic ray events by the presence of two energy deposits above 30 MeV in 
the outermost calorimeter layers [20]. A fraction about 80% of the cosmic ray events 
are identified and rejected at the trigger level with this technique. Further suppres- 
sion of the DA<I>NE background events and cosmic rays is performed by an off - line 
filter called FILFO (FILtro FOndo: background filter). FILFO identifies background 
events at a very early stage of the data reconstruction using only information from the 
calorimeter [47]. 

For the search of the Ks 37r° decay only the calorimeter signals are used to trigger 
the event. Two energy deposits above threshold about 50 MeV for the barrel and about 
150 MeV for the endcaps are required [16]. 



^As it was mentioned a part of the e+e ^ e+e events are gathered for detector monitoring and 
caUbration. 
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Chapter 4 

First stage of the event selection 



The meson produced in the e~^e~ collision at DA$NE is in a pure J^*" = 1 state. 
Since the — > KsKl decay is driven by the strong interaction, the initial KsKl state 
is antisymmetric with the same quantum numbers and can be written in the rest frame 
as: 

\i) = N- [\Ks{p)) \KL{-fP)) - \Ks{-m , (4.1) 

where p denotes the momentum of each kaon and is a normalization factor [35]. 
Since the resonance is moving with a small momentum in the horizontal plane ~ 
13 MeV/c Ks and mesons are produced almost back - to - back in the laboratory 
frame. Therefore, observation of a (Ks ) decay ensures the presence of a Ks (K^ ) 
meson travelling in the opposite direction^ [48]. Thus, at DA$NE we obtain pure 
Ks and Kl „beams" with precisely known momenta and flux, which can be used to 
measure absolute branching ratios [20]. In this chapter the Ks tagging technique with 
the detection of the K^ interaction in the KLOE calorimeter is described. 

4.1 Identification of Ks via detection of Ki 

Neutral kaons produced at KLOE have a velocity in the rest frame equal to (3 ^ 0.22. 
This corresponds to the K^ time of flight from the interaction point to the calorime- 
ter equal to about 31 ns, which means that about 60% of produced K^ mesons reach 
the calorimeter without decaying [20]. Kl mesons interact in the calorimeter with an 
energy release up to ~ 497 MeV (so called „Kl - crash"). Thanks to the exceptional 
timing capabilities of the KLOE calorimeter^ and the low velocity of kaons one can 
use the Time of Flight technique to tag the Ks meson, as described in the next section. 

' We refer to the process of defining a Ks or Kl sample as tagging; observation of a Kl {Ks ) 
decay tags the presence of a Ks (Kl ) meson and allows for the determination of its momentum [48]. 

^ For an energy release of 100 MeV the resolution of time measured by the calorimeter amounts to 
about 0.3 ns, which corresponds to about 1% accuracy in the determination of the Kl velocity [20]. 
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Adding the information about the position of the energy release (K^ cluster), the di- 
rection of the Kl flight path can be determined with ~ 1° angular accuracy [20]. This 
allows to estimate the Kl momentum vector and as a consequence, knowing the four 
- momentum, to determine the four - momentum of the tagged Ks meson. 

4.1.1 Identification of the Kl meson 

The identification of the Kl interaction in the calorimeter is performed after tracks 
reconstruction and association to the clusters^, and after the preselection aiming at the 
rejection of events with Kl meson decay inside the drift chamber. Events for which 
there is one reconstructed vertex with two tracks having opposite curvature are rejected. 
Moreover, an event is discarded if there are two reconstructed tracks, having opposite 
curvature, associated to two vertices reconstructed less than 30 cm away from the in- 
teraction point in the transverse plane. These cuts reject most of the background events 
with Kl decaying before reaching the calorimeter [49] . For each surviving event we 
look for the Kl clusters in the calorimeter taking into account only clusters not associ- 
ated to any track. For each that kind of clusters we calculate velocity of the contributing 
particle defined in the laboratory frame as: 



where Rd denotes the distance from the e"^e~ interaction point to the reconstructed 
position of the cluster center, td stands for the measured time of flight of the particle 
and c is the speed of light. Since in the (f) rest frame kaons have a well known velocity 
/3 ~ 0.22 it is convenient to transform (3ci to this reference frame: 



where /^^ denotes the velocity of the meson in the laboratory frame and a stands for 
the angle between the momentum vector and a direction vector connecting the inter- 
action point with the cluster position. The distribution of (3cr is presented in Fig. 4.1a. 
A big peak around /3cr = 1 corresponds mainly to clusters formed by gamma quanta 
from Ks mesons which decay very close to the interaction point. Clusters originating 
from the decay products of remaining Kl decays are instead characterized by smaller 
velocities distributed in the range from (3cr ~ 0.28 to (3cr ~ 1 [49]. In Fig. 4.1a one can 
also see a structure around (3cr ~ 0.22 corresponding to Kl mesons and a smaller peak 
for very low velocities. The latter peak originates mainly from the charged particles 

^ The track - to - cluster association procedure establishes correspondence between tracks in the 
drift chamber and clusters in the calorimeter. 




(4.3) 



cr 



1 + (3<pPdCosa 



4.1 Identification of Ks via detection of Kl 
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Figure 4.1: Distribution of the (3cr velocity reconstructed for clusters not associated 
to the track for a sample of events before identification of the Kl interaction in the 
calorimeter. The spectra are made for all clusters before (a) and after (b) cut on energy 
E,r > 100 MeV. 



(e.g. pions) for which the track - to - cluster association procedure failed. Additionally, 
most of these clusters are characterized by energy deposits smaller than 100 MeV [49]. 
The (3cr distribution after the cut on Ecr > 100 MeV is shown in Fig. 4.1b, where 
one can see a big suppression of the clusters with lowest velocity. Therefore, clusters 
originating from the interaction in the calorimeter are defined with the following 
conditions: 

0.17 < I3cr < 0.28 

Ecr > 100 MeV , (4.4) 
where Ecr is the energy of the cluster. 

The main remaining background sources to this tagging algorithm are the cosmic 
muons entering KLOE through the intersection between the barrel and endcap calorime- 
ters. Such muons may give a signal in the calorimeter without a track in the KLOE drift 
chamber. The other contributions to the background originate from DA<I>NE activity 
and — )■ K^K~ decays [43]. The angular momentum of the K^Ks system is equal to 
the spin of the meson s = 1. Therefore, kaons from the decay are mostly emitted 
in the direction perpendicular to the beam axis and the background can be additionally 
suppressed selecting only „Kl - crash" clusters in the barrel [43]. 
Kl meson interacting in the calorimeter usually induces more than one cluster, there- 
fore to estimate the direction of the tagging Kl meson we consider the „fastest" can- 
didate cluster which was produced as the first one. 
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4.1.2 Kl - momentum estimate 

In the (j) — )■ KsKi decay the four - momentum can be determined completely 
knowing the center of mass energy y/s, the cj) momentum vector and a angle be- 
tween the momentum and the Kl flight direction determined in the laboratory frame 
from the reconstructed center of the cluster. For ^/s and we use the mean values 
measured for each running period using the gathered sample of e"'"e~ scattered at large 
angles. Determination of Pk^^ allows to calculate the four - momentum of the tagged 
Ks meson: Pr, = P<^ - Pr,. 



Chapter 5 

Normalization sample 



Registration of the Kl interactions in the calorimeter allows for the simultaneous iden- 
tification of the Ks meson. Since one of the goals of this measurement is to determine 
the Ks 37r° decay branching ratio, the number of events selected as the signal has to 
be normalized to the number of all Ks decays. To this end the Ks — ?• 27r° — )• 47 events 
were also counted (further on they will be referred to as the normalization sample). 
This process is one of the main Ks decay channel with well - known branching ratio 
BR{Ks 27r°) = 0.3069 ± 0.0005 [19]. The number of events produced for both the 
signal and the normalization sample can be expressed as: 

= L-a^- BR{4> ^ KsKl) ■ BR{Ks ^ 27r°) ■ ■ e,, (5.1) 
iVs. = L-a^- BR{(t) ^ KsKl) ■ BR{Ks ^ 37r°) ■ £3. ■ e,, , 

where L is the integrated luminosity, denotes the total cross section for produc- 
tion, stands for the tagging efficiency and e-i-,, and e2T, are the identification efficien- 
cies for the appropriate channel. The ratio: 

^ BRjKs ^ 37r°) ■ 63. 
iV2. BR{Ks ^ 27r0) -62^ 

allows for the BR{Ks — )■ 37r°) determination independently of L, a^, ecr and the 
(p — )■ KsKl branching fraction avoiding all the systematic effects originating from 
measurements of these quantities. 

After identification of the K^ meson interacting in the calorimeter the preselection is 
based on the number of reconstructed 7 quanta in each event. To this end we consider 
only calorimeter clusters not associated to any track reconstructed in the drift chamber. 
Moreover, the reconstructed time of the cluster t^i should be compatible with the time 
of flight of photon equal to Rd/c, where R^ is the distance from the cluster position 
to the interaction point' and c denotes the speed of light. We assume the two times to 

' Since the mean free path of the Ks originating from the 4> decay amounts to about 6 mm (the kaon 
velocity /? ^ 0.215) , which corresponds to negligible time delay equal to ^ 100 ps, we assume that it 
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Figure 5.1: The experimental distribution of the reconstructed 7 quanta mul- 
tiplicities after imposing the tight A'5 tag requirements ( E^r > 150 MeV and 
0.200 < Per < 0.225 ) and acceptance cuts defined by Eqs. 5.3, 5.5 and 5.6. 



be consistent if: 

\tci - Rci/c\ < MIN(3.5 ■ at{E,i), 2 ns) , (5.3) 

where at is the calorimeter time resolution parametrized as a function of the cluster 
energy E^i: 

ME^i) = , © 140 ps. (5.4) 

The cutoff on 2 ns is used to reduce the number of the machine background clusters 
accidentally overlapping with the event. To this end we apply also cuts on the minimal 
cluster energy and polar angle: 

E,i>7MeV (5.5) 



I cos{9ci)\ < 0.915 23.8° < O^i < 156.2° . (5.6) 

Distribution of the 7 quanta multiplicities is shown in Fig. 5. 1 . At this stage of analysis 
we select two data subsets: the signal sample which consists of events with six recon- 
structed photons and the normalization sample of Ks 27r° candidates with N-y = 4. 
For both channels the expected background as well as the detector acceptance and 
the analysis efficiency is estimated using the Monte Carlo simulations based on the 
GEANT3 package [50]. The simulations incorporate a detailed geometry and mate- 
rial composition of the KLOE apparatus and all the conditions of the experiment e.g. 



decays exactly in the interaction point. 



5.1 Comparison between data and simulations results for 47 events 
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DA$NE background rates, position of the interaction point and beam parameters^. 

5.1 Comparison between data and simulations results 
for 47 events 

For the search of rare processes like the Ks — 37r° the estimation of the background 
has to be as precise as possible. In our research to this end we use Monte Carlo simula- 
tions described briefly in the introduction to this chapter. Moreover, the determination 
of the efficiencies of cuts and discriminant analysis is also based on the simulated sam- 
ples of events. Therefore, we have checked the reliability of the KLOE Monte Carlo 
simulations and optimized them for the best possible description of the experimental 
data. 

Since the reconstruction efficiency of clusters in the calorimeter is slightly higher for 
simulations compared to the measured data we apply a correction determined based 
on the —7- 7r+7r^7r° sample^. The efficiency for both data and simulations were 
parametrized as a function of the 7 quanta energy and polar angle. The correction 
is then applied deleting randomly photons from the simulated events with a probability 
equal to the ratio of efficiencies for data and simulations [51, 52]. 
Apart from the cluster reconstruction efficiency the simulations were corrected also 
for the energy scale of the reconstructed gamma quanta. The necessity of this addi- 
tional correction is justified in Fig. 5.2a where we observe a small shift between the 
distributions of the reconstructed Ks mass for data and simulations before the correc- 
tion. The procedure of the energy scale correction for Monte Carlo events is based 
on the Ks — )• 27r° — )• 47 sample which is almost background free. For both, data and 
simulations the following variable has been constructed: 

^E, = , (5.7) 

Efit 

where Ej is the energy of reconstructed gamma quantum and Efn denotes energy 
of the same gamma quantum corrected by the kinematical fit procedure. The gamma 
quanta were then divided into groups of 20 MeV with respect to the Efu. For each 
group the AE^ distribution was fitted with the Gauss function. The mean values of the 
fitted Gauss distributions are shown in Fig. 5.3a. As it can be seen Monte Carlo simu- 
lations systematically underestimate the data. The energy scale correction was applied 
by increasing the cluster energies by a factor parametrized as a function of Efu. For 

^The detailed description of the KLOE Monte Cario simulation program GEANFI can be found in 
Ref. [44]. 

■'The detailed description of the cluster reconstruction efficiency studies for data and simulations can 
be found in Ref. [52]. 
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Figure 5.2: Comparison between the reconstructed Ks mass distributions for data (red 
points) and simulations (blue histogram) before a) and after b) the energy scale correc- 
tion. The Ks mass is reconstructed from the Kg — )■ 27r° events. 
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Figure 5.3: Distributions of the mean A^E^ for data (black squares) and simulations 
(blue circles) without energy scale correction a) and after the correction b). The values 
are obtained using the fit described in the text. 
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the first 20 MeV bin we started with 2.4% shift while for every next group of clusters 
the correction was decreasing by a factor of 0.1 %. The result is presented in Figs. 5.2b 
and 5.3b where one can see much better agreement with the data. 
For further validation and tuning of the simulations, after applying the cluster effi- 
ciency correction, we have determined relative fractions of number of events with 
given 7 quanta multiplicity k with respect to the total number of events with 3-6 
reconstructed photons: 

F.^^^, (5.8) 

and compared experimental values of F^. with results of simulations. The distributions 
of the relative fractions as a function of the KLOE running period are presented in 
Fig. 5.4. The agreement between data and simulations is reasonable apart from the 
most important multiplicity Fq where results of the simulations systematically over- 
estimate the data during the whole data taking period. To understand the source of 
this discrepancy we have determined the probabilities to find one ( P^i ) or two ( Pa2 ) 
accidental clusters in the prompt time window defined in Eq. 5.3 for both data and sim- 
ulations. To this end we have considered clusters in so called early time window i.e. 
background clusters originating from earlier bunch crossing fulfilling the condition: 

(tez - Rci/c) e [-54, -14] ns , (5.9) 

which corresponds to about fifteen groups of accidental clusters sources'^. The times of 
these clusters were then shifted by a number of bunch crossing periods T^/ obtained 
for each event from the time of the earliest cluster t^^ as follows: 

= U + T,^ ■ INT(^^;^) , (5.10) 

and next the acceptance cuts defined in Eq. 5.3 were imposed. This allowed to deter- 
mine the fraction of events with one or more accidental clusters in the acceptance and, 
as a consequence, to calculate P41 and P42 (see Tab. 5.1). 

The values of P41 and Pa2 were next used to estimate the probability for a cluster to 
produce one ( Ps\ ) or more fragments ( P52 ) reconstructed as an additional cluster. 
In order to do that the true relative photon multiplicities F^''"^ not affected by the ac- 
cidental activity or cluster splitting^ were determined based on the information about 
the decay chain of the simulated events. Knowing the F^^^ values and the determined 

^The minimum bunch crossing period of DA$NE is equal to T, / = 2.715 ns. 
^ In the simulations the accidental clusters can be ignored referring to the GEANT KINE indices for 
particles contributing to each cluster 
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Figure 5.4: Comparison of the distributions for data (black squares) and background 
simulations (blue circles). The solid and dashed lines denote linear fits (F^ = const.) 
to the spectra of simulations and data, respectively. The results of the fit are reported 
on the top of each distribution. 
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Pai [%] 


Pa2 m 


Psi [%] 


Ps2 m 


DATA 


0.378 ± 0.004 


0.025 ± 0.001 


0.30 ± 0.01 


0.0103 ± 0.0001 


SIMULATIONS 


0.492 ± 0.004 


0.027 ± 0.001 


0.31 ±0.01 


0.0156 ± 0.0002 



Table 5.1: The probabilities to find one ( Pai ) or two ( P42 ) accidental clusters and 
to reconstruct one ( Psi ) or more ( P52 ) splitted clusters estimated using events in 
out-of-time window and fit to the distributions, as it is described in the text. 



probabilities of accidental coincidence P^i and P42 we can fit the measured dis- 
tributions treating the Psi and Ps2 as the free unknown parameters^. The obtained 
splitting probabilities are presented for both data and simulations in Tab. 5.1. These 
results show that for the simulations there is about 50% more events with two split- 
ted clusters and about 30% more events with one cluster originating from the machine 
background, which explains the discrepancy for the 6 - gamma events. The technique 
used to account for this difference is presented in chapter 6. 

As described in chapter 4 at the preselection stage we cut on the velocity (3cr and ener- 
gy Ecr of the Kl meson, therefore the simulations of its interaction in the calorimeter 
should be also realistic and precise. The comparison of the /3cr and Ecr distributions 
for data and Monte Carlo simulations is presented in Fig. 5.5a and 5.5b. It is clear, 
that the simulated Kl velocity is in a reasonable agreement with data while there is a 
big discrepancy in the Ecr distribution. Therefore before the cuts optimization another 
correction for the Monte Carlo simulated events had to be applied. A small correction 
to the (3cr was made adding a shift of 1 % multiplied by a Normal - distributed random 
number. The Ecr was instead modified by 2.5% correction increasing every 1 MeV by 
0.03%. The resulting distributions are shown in Fig. 5.5c and 5.5d where one can see 
a much better agreement with data. 

Finally applying all the corrections described before we have compared some other 
simulated inclusive distributions for the 4 - gamma sample with the experimental ones. 
As it can be seen in Fig. 5.6 the agreement is reasonable so we can proceed with further 
analysis and counting of the Ks 27r° events. 



''The detailed description of the fit procedure and used probabilistic model is described in Ref. [52], 
where the technique to measure the Pai and Pa2 probabilities are also presented more detailed. 
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Figure 5.5: Distributions of energy and velocity of the tagging meson for data (red 
points) and simulations (blue histogram) before ( a and b ) and after the corrections ( c 
and d ). Here only events with four gamma quanta are taken into account. 
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Figure 5.6: Inclusive distributions for the Ks 2tx^ — > 47 decay after all corrections 
applied to the Monte Carlo simulations (blue histogram). The red points denote the 
experimental result. 
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5.2 Background estimation and counting of the Ks 
27r° — )► 47 events 

Before the final counting of the Ks — )■ 27r° sample one has to estimate the number of 
events originating from processes different than Ks 2%^ for which we have found 
four reconstructed gamma quanta. To this end we have used the Monte Carlo simulated 
events of (p decays for which the true decay chain is known. Based on this information 
we have found that the background events constitute a negligible fraction of the 4 - 
gamma sample amounting to about 0.1%^. This allows for counting of the Ks — )■ 
27r° candidates without any further cuts. With the tight requirements for reconstructed 
Kl energy ( E^r > 150 MeV and velocity (0.200 < (3„. < 0.225 ) we have found the 
following number of events with four reconstructed gamma quanta^: 



2-K 



(7.533 ±0.018) X 10^ . 



(5.11) 



The Monte Carlo simulations allow also for the determination of the Ks — 27r° — 47 
selection efficiency. To this end we have used a simulated sample of Ntot = 59141 
Ks — )■ 27r° events fulfilling above mentioned tight conditions for f3cr and Ecr- Next 
we apply the acceptance cuts and count events with a given number of reconstructed 
gamma quanta (see Tab. 5.2). The efficiency is then defined as: 



Ajrec 



(5.12) 



tot 



T\jrec 



denotes the number of events with four reconstructed gamma quanta and Ntot 
stands for the total number of simulated Ks — t- 27r° events. In both cases we count 
only events fulfilling the tight Kl tag conditions. Taking into account the estimated 



k 





1 


2 


3 


4 


5 


6 


7 


8 




25 


181 


3059 


16185 


39012 


636 


40 


1 


2 




0.043 


0.31 


5.17 


27.37 


65.97 


1.08 


0.068 


0.0017 


0.0034 



Table 5.2: The number of events Nl^'^ with a given multiplicity of reconstructed clus- 
ters k for Ntot = 59141 simulated KsKl 4:jKl events with 0.200 < /3cr < 
0.225 andEc,. > 150 MeV. 

selection efficiency: e2,r = 0.660 ± 0.002stat we can determine the final number of 
produced events: 



N2^ 



;i. 142 ±0.005) ■ 10^ 



(5.13) 



This number will be used for the normalization of the Ks — )■ 37r° branching ratio. 



^These are above all the </> — > K^K events. 

^ The justification of using this set of cuts will be presented in chapter 6. 
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Chapter 6 

Search for the Kg Stt^ 67 signal 



After the Ks tagging via the Kl interactions in the calorimeter the selection of the 
Ks 37r° decay is performed by searching for six photons from the decay of the 
pions. After the preselection with the conservative^ Kg tag requirements 
( Ecr > 100 MeV and 0.17 < /3cr < 0.28 ) we have found 76689 events with six re- 
constructed 7 quanta. For these events we perform further discriminant analysis to 
increase the signal to background ratio. As it was in the case of the Ks — 27r° normal- 
ization sample we have considered only calorimeter clusters not associated to any track 
reconstructed in the drift chamber and imposed acceptance cuts defined with Eqs. 5.3, 
5.5 and 5.6. 

As it was mentioned in chapter 5, the expected number of background events as well 
as the analysis efficiency is estimated using the KLOE Monte Carlo simulations. All 
the processes contributing to the background were simulated with statistics two times 
larger than the measured data. Moreover, for the acceptance and the analysis efficiency 
studies the dedicated Ks — )■ Stt" signal simulations were performed. The signal is 
generated taking into account the conditions of the KLOE experiment and assuming 
branching ratio equal to the best known upper limit [16] increased by a factor of 30. 
Background to the searched Ks 37r° — 67 signal originates predominantly from 
the Ks 27r'^ events. The four photons from this decay can be reconstructed as six 
due to fragmentation of the electromagnetic showers (so called splitting). These events 
are characterized by two clusters with low energy and position near the place where 
the true photon hits the calorimeter. The additional fake clusters can be generated also 
by the accidental coincidence between the (p decay event and the DA$NE background. 
The other source of background are (p — )■ KsKl — )• 7r+7r^,37r° events. Methods of 
suppressing this kind of background will be discussed in the next section. 



^ Such requirements have been used in the previous analysis of the KLOE data [16]. 
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6.1 Discriminating variables and the signal region def- 
inition 

6.1.1 Rejection of the — > Stt^ events 

The most dangerous source of background for our analysis are the — )■ KsK^ 
{Ks — > vr+vr", — Stt^) events. In this case one of the charged pions can interact in 
the DA$NE low - beta insertion quadrupoles producing neutral particles which may 
ultimately simulate the signal of Kl interaction in the calorimeter, while the Kl meson 
decays close to the interaction point and produces six photons. To suppress this kind 
of background we first reject events with charged particles coming from the vicinity 
of the interaction region. Since Ks decays near the interaction point we reject events 
with at least one track satisfying: 

Pt = \l x^PCA + VpcA < 4 cm (6.1) 
\zpca\ < 10 cm , (6.2) 

where PCA denotes the point of the closest approach of the reconstructed trajectory 
to the interaction region. In principle one could reject all events with the track recon- 
structed in the drift chamber. However, this would decrease the reconstruction effi- 
ciency of the Ks — 37r° signal because of the relatively high probability of an ac- 
cidental coincidence of the real event with tracks generated by DA$NE background. 
The choice of conditions defined in Eq. 6.1 minimizes the loss of the signal events due 
to this effect. 




Figure 6.1: Simulated distributions of the reconstructed K^ energy E^r and velocity 
[5cr for the (p KsKl {Ks Stt", Kl) 67 signal (solid) and -)■ KsKl 
{Ks — TT+TT^, Kl — Svr^) 7r"'"7r~67 background events (dashed) after the conser- 
vative cuts on reconstructed (3cr and E^r described in Eq. 4.4. 
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Figure 6.2: The rejection efficiency for the — KgKi — t- tt+tt", Stt^ reaction chain 
after cutting events with tracks (TrV), after the tight requirements for Kl interaction 
in the calorimeter (HardKI) and both cuts (TrVH). The result was obtained based on 
the simulated background sample. 

As it is presented in Fig. 6. 1 the distribution of the velocity /5cr of Kl candidates recon- 
structed from signals in the calorimeter for events with fake „Ki - crash" is relatively 
flat, while there is a clear peak around 0.215 for the signal events. Similarly there is 
a substantial difference in the reconstructed energy spectrum of the interacting Kl 
mesons and fake „Kl - crash" events. Thus, tightening the cuts on E^r and f3cr allows 
to reject almost all events belonging to this background category: 



These cuts reduce the Ks tagging efficiency from 34% to 23% but we gain a big re- 
duction factor (about 60) on the most important source of background (see Fig. 6.2). 

6.1.2 The kinematical fit 

In the next stage of the analysis we select only kinematically well defined events. To 
this end we perform the kinematical fit procedure based on the least squares method 
with the following set of variables as an input: 

• The total energy of the system ^/s 

• The momentum and decay vertex of the meson 

• The four - momentum vector and decay vertex of Ks meson 



> 150 MeV 



0.200 < I3^r < 0.225 . 



(6.3) 
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Figure 6.3: a) distribution of the kinematical fit for 47 events of data (red pionts) and 
simulations (blue histogram); b) The distributions from kinematical fit for simulated 
Kg —7- 37r° signal events (solid) and background (dashed). 



• Energy of each 7 quantum together with its time and position measured in the 
calorimeter 

The initial value of the Ks four - momentum vector is determined for each event using 
the position of the reconstructed cluster while for ^/s and we use the mean 
values measured for each running period using the gathered sample of e^e~ scattered 
at large angles. Using this set of variables we can construct the following quadratic 
form: 

= (Y- YofV~\Y - Yo) + A ■ G(Y) , (6.4) 

where Yq is a vector of measured variables, V denotes covariance matrix, Y denotes 
the vector of corrected variables fulfilling the kinematical constraints G(Y) and A is a 
vector of the Lagrange multipliers. The quadratic form is then minimized with respect 
to Y and A leading to the determination of the best corrected variable values. The min- 
imum value of can be treated as the - like variable with a probability function 
that can be used for hypothesis testing and evaluation of the goodness of the fit. As 
the constraints we impose the total four - momentum vector conservation, Ks mass 
requirement and consistency of the time of flight determined from the cluster position 
of each 7 quantum with its time reconstructed in the calorimeter. The developed algo- 
rithm was first applied to the kinematically well defined Ks 27r° 47 events to 
test how reliable is the fit procedure. The distributions for data and simulations for this 
kind of events are shown in Fig. 6.3a. The same distributions for generated background 
events and Ks 37r° signal are shown in Fig. 6.3b. It can be seen that applying a cut 
on around 40 we obtain a small background rejection, about 30%, with a good sig- 
nal efficiency (~70%). This does not improve significantly the signal to background 



6.1 Discriminating variables and the signal region definition 



41 



ratio but allows to reject bad quality reconstructed events. 



6.1.3 Testing the 37r^ and Itt^ hypotheses 

In order to reject events with split and accidental clusters we look at the correlation 
between two - like discriminating variables xL ^i^^ xlw- xL calculated by an 
algorithm selecting four out of six clusters best satisfying the kinematic constraints 
of the two-body decay, therefore it verifies the Ks — )■ 27r° — 47 hypothesis. The 
pairing of clusters is based on the difference between reconstructed 7r° masses M^o 
and M^o with respect to the PDG value Mpdg [19] and on the opening angle of the 
reconstructed pions trajectories in the Ks center of mass frame 9-„t, which should be 
equal to 180° for the Ks — J- 27r° events. Moreover, we check the consistency of the 
determination of the Ks four - momentum vector P^^^ . It is performed by comparing 
the Pj^-g determined from the reconstructed four - momentum of K^with the sum of 
the 7 quanta four - momenta Prec = Ylt=i ^7i- every possible pairing choice the 
algorithm calculates the xL defined as: 



(M,o - MpDG? (M,o - MpDG? {6^^ - Tlf 



~ „2 ~^ 



2 ' _2 ' _-2 ' _2 

{pks - E p^) ' [pk - E pi) ' {pks - E p^i ' 

CTp CTp Up 

-i X ^ y -« 2 



(6.5) 



The minimization of the xL gives the two photon pairs which out of the six clusters 
fulfills best the criteria expected for the Ks — )• 27r° — j- 47 hypothesis. The resolutions 
used in Eq. 6.5 were estimated using the Ks 27r° — )• 47 normalization sample. The 
four - momentum vectors of the reconstructed pions were used to make distributions of 
the differences used for the xL calculation (see Fig. 6.4), which were then fitted with 
the Gauss functions allowing for the determination of the variances used in Eq. 6.5. 
All the values of the parameters are gathered in Tab. 6.1. 

The Xstt instead verifies the signal hypothesis by looking on the reconstructed masses 
of three pions. For every choice of cluster pairs we calculate the quadratic sum of 
the residuals between the nominal 7r° mass and the invariant masses of three photon 
pairs [21]: 

xL = — 2 + — 2 + — 2 • (6-6) 

(^L ^l. 

As the best combination of cluster pairs we take the configuration with lowest Xa^- 
The resolution of pion mass a^T, was estimated applying the algorithm to the simulated 
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PARAMETER 


DATA 


SIMULATIONS 




18.66 ±0.04 


19.24 ±0.04 




18.84 ±0.04 


19.36 ±0.04 




44.07 ± 0.08 


46.83 ± 0.09 




25.93 ± 0.05 


27.55 ± 0.06 




26.12 ±0.05 


27.55 ± 0.06 




23.48 ± 0.05 


24.39 ± 0.06 




0.1238 ±0.0002 


0.1257 ±0.0002 




17.0 ±0.5 


17.0 ±0.5 



Table 6.1: List of parameters used in the calculation of xItt- The values were 

obtained by fitting Gaussian functions to the distributions presented in Fig. 6.4 and to 
the distributions of pion masses reconstructed by minimization of the value using 
the simulated signal events. 




Figure 6.4: Distributions of parameters used in calculation for data (red points) 
and the simulations of the Ks Itx^ channel (blue histogram) after the energy scale 
correction described in Sec. 6.2. 
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signal events (see Tab. 6.1). In the definition of ^° ^^^e into account the 

difference between Fk^ and Free = Ylt=i ^^tj because in this case it is the same for 
each combination of photon pairs. 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 



%3jt ^CSti 

Figure 6.5: Scatter plot of xlir versus xL simulated Ks 37r° signal (left) and for 
the background after rejection of the — )■ KsKl vr"^7r~, 37r° events (right). 

As it is presented in Fig 6.5 the Ks — 37r° signal is characterized by low values of 
Xstt and relatively high values of xlw Background events are instead spread on the 
large area of (xL? xln) plane with a maximum at low xlw being well distinguishable 
from the signal. Nevertheless in the region populated by signal we find also some 
background events, mainly the (J) — KsK^ — )■ 7r+7r~,37r° category even though it 
is already strongly suppressed by cuts defined with Eqs. 6.1 and 6.3. Further analysis 
dedicated to the rejection of the background due to the Kg — )■ 27r° decay is described 
in the next section. 

6.1.4 Improvement of the Ks — > 27r° background suppression 

Since we are looking for a very rare decay and have to deal with a large background 
the rejection power of the data analysis has to be as high as possible. The two addi- 
tional four - momentum vectors of photon candidates reconstructed based on clusters 
originating from the machine background or shower fragments results sometimes in an 
invariant mass close to the mass of vr^. Thus, the Xstt algorithm does not distinguish all 
the Ks — 27r° decays from the signal and we need another discriminant variable to 
refine the background rejection. 

Events with two accidental clusters can be identified by measurement of the difference 
Ai? between the Ks energy determined from the reconstructed K^ four - momentum 
and the sum of energies of four gamma quanta selected by the xL algorithm. For the 
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Figure 6.6: Distributions of variables used to refine the rejection of the Ks — 2n^ 
events. Dashed histogram indicates simulated background events from Ks — ^ 27r° and 
solid histogram represents the Monte Carlo simulations of the signal. The distributions 
are made for events with < 300. The variables are described in the text. 



Ks — !■ 27r° background this variable is close to zero since the event is kinematically 
closed. On the other hand for the Ks — )■ 37r° events AE should peak around 135 MeV 
since the rejected clusters have an energy comparable to the pion mass. In order to 
make the cuts as independent as possible of the energy resolution determination we 
use the normalized AE: 

AE/aE = , (6.7) 

where the used value of aE is equal to the one listed in Tab. 6.1. An example of the 
AE/aE distributions for simulated background and signal are presented in the left 
panel of Fig. 6.6. One can see that cutting around 1.8 allows to reject about 60% of 
background events keeping the signal efficiency at the level of around 80%. 
Further on, events with splitted clusters are suppressed with cut on the distance be- 
tween center of reconstructed clusters. Here we take advantage of the fact that the dis- 
tance between splitted clusters is on average smaller than the distance between clusters 
originating from 7 quanta of the Ks — )■ 37r° decay. For all possible pairs of clusters 
(z, j) we calculate the distance: 

Rij = \l {xi - XjY + [vi - VjY + {zi - ZjY , (6.8) 

where (xj, yi, Zj) and (xj, yj, zj) are the position coordinates of clusters reconstructed 
in the calorimeter. We look then for the minimum of obtained values Rmin = MIN {Rij ) 
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and impose the following cut: 

Rmin > 65 cm . (6.9) 

This suppresses the background by about 30% retaining about 85% of signal events 
(see right panel of Fig. 6.6). The AE/as and Rmin close the cuts sequence improving 
the background rejection power of the analysis chain. 

6.2 Background estimation 

Search for the rare decays like Ks — Sn^ demand a precise knowledge of the back- 
ground processes. To this end we used realistic Monte Carlo simulations based on 
GEANT3 package. As it was mentioned in chapter 5 there is a substantial discrepancy 
between the data and results of the simulation for the six photons sample, which had to 
be fixed. To this end we first divide the background simulations into categories using 
for each event information about the decay chain and about particles contributing to 
each cluster. All the categories are next used in the fitting procedure to the experimental 
data, which allowed for the determination of the weighting factors for events belong- 
ing to each class of simulated background. The procedure used to determine the event 
weights is described in the following subsection, where we present also small correc- 
tions related to the simulated Rmin distribution and rejection of events with charged 
particles. 

6.2.1 Event weights determination 

A below listed classes of background events were recognized^: 

• Fakes: the — )■ KsKl {Ks vr+vr", Kl 37r°) events together with non 
- KsKl channels like K+K~ or -)> tt+tt^tt^ (about 2%) 

• 2A-I-1A1S: events with two accidental clusters (about 24%) or events with one 
accidental and one splitted cluster (about 6%) 

• 2S: events with two splitted clusters (about 60%) or with more than two acci- 
dental or splitted clusters, as well as Ks 2n^ events with only one splitted or 
only one accidental cluster (about 8%) 

The category called Fakes contains, apart from the standard — KsKl — {Ks 
7i^ti~,Kl 37r°) events, small admixture of the non - KsKl background, which 
is however suppressed by the ^^id AE/aE cuts. As regards other classes of back- 
ground it is relatively easy to justify the presence of more than two clusters originating 

^ The given fractions of events refer to the simulated sample after the preselection described at the 
beginning of this chapter 
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Figure 6.7: The X2n versus xI-k scatter plots for different background categories. The 
label of each plot is defined in the text. The corresponding experimental spectrum is 
also shown. 



from the machine accidental activity or from cluster splitting. 

The Ks — 27r° events with only one splitted or only one accidental cluster in princi- 
ple should not pass the requirement of the registration of six gamma quanta. However, 
more detailed analysis revealed, that the additional real cluster can be generated by the 
initial state radiation or by the wrongly reconstructed cluster^. For both data and 
all the background categories we have made the (xItt'XL) scatter plots. 
The Monte Carlo distributions were then fitted to the data as a linear combination pro- 



It has turned out, that for many events in this category the additional cluster is very close to the place 
where Kl interacted with calorimeter Thus they can result from wrongly reconstructed Kl showers. 
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viding scaling factors Wt for each background category used next to weight events: 

Data(xL, Xl) = ■ Sim^^(xL, xl) + W^^^^'^^' ■ Sim^^+^^^S(xL, xl) 
+ iy^^'^--Sim^^'^-(xL,xL), 

where Data(x27r; xlw) denotes the experimental (xL; xln) scatter plot and 

Sim(x|^, xItt) stands for the simulated distribution of each background category. The 



Category 


Nmc 


Nfit 


Wt 


2A+1A1S 
2S 
Fakes 


64544 ± 254 
145996 ± 382 
5670 ± 75 


26346 ± 236 
43446 ± 283 
6897 ± 148 


0.4082 ± 0.0040 
0.2976 ±0.0021 
1.216 ±0.031 



Table 6.2: Scaling factors for Monte Carlo background categories used in the fit to the 
data. Nmc and Nju denote the number of events in each category before and after the 
fit, respectively. 

results of the fit are gathered in Tab. 6.2. The quality of the procedure used to refine the 
simulations for the six - photon sample can be controlled by the comparison of simu- 
lated and experimental inclusive distributions of discriminating variables described in 
previous sections just after preselection. This comparison is shown in Fig. 6.8 and im- 
plies that the agreement between data and background simulations is reasonable after 
the fit. Another check was done dividing the (xL; xI-k) plane onto different regions 
to compare the experimentally observed number of events with expectations based on 
the Monte Carlo simulations. Five control boxes were chosen around the signal region 
defined with preliminary cuts on xlw and xlw (see Fig 6.9). The number of events reg- 
istered in each box is reported in Tab. 6.3, where one can see that simulations results 
agree with experiment when taking into account the statistical uncertainties calculated 
using the standard deviations of the scaling factors (see Tab 6.2): 

ANho:c = \J{Wfy/N^Y + (iy2A+lAlS^^2A+lAlS)2 + Fakes ^^Fakes^ ^ 

(6.10) 

where denotes the number of events belonging to each background category. 
As it can be seen in Fig. 6.10a, despite of application of the weighting factors for 
different background categories there is still a small discrepancy between data and 
simulations for the Rmin discriminant variable. Moreover it turns out, that we reject 
more events with charged particles having f3cr > 0.220 (see Fig. 6.10c). Therefore, 
we introduce a small correction shifting the Rmin by 2 cm smeared with a Normal- 
distributed random number. The Rmin distribution after this correction is presented 
in Fig. 6.10b where we observe much better agreement with data. For the simulated 
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Figure 6.8: Inclusive distributions of discriminating variables for data (red points) and 
the simulations of the background (blue histograms) weighted with Wt factors as de- 
scribed in the text. 





SBOX 


DOWN 


UP 


CDOWN 


CUP 


CSBOX 


DATA 


200 ± 14 


416 ±21 


7 ± 3 


14385 ± 120 


16321 ± 128 


17634 ± 133 


MC 


228 ± 10 


313 ± 12 


8±3 


14380 ± 134 


16143 ± 124 


17940 ± 123 



Table 6.3: The number of events populating control boxes in the {xln^ xL) plane de- 
fined in Fig. 6.9 after tight requirements for the reconstructed Kl energy and velocity. 
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Figure 6.9: Zoom of the (xItt; xL) distribution for data illustrating the division of the 
plane into control boxes. 

events with /3cr- > 0.220 belonging to the Fake category a correction based on the „Hit 
& Miss" method was applied. For each rejected simulated event we draw a random 
uniformly-distributed number from the range (0.; 1.), if this number was greater than 
0.93 the event was passed for the further analysis. The /3cr distribution for the rejected 
events after the correction is presented in Fig. 6.10d. 



6.3 Optimization of selection criteria 

After all corrections described in last section as well as in chapter 5, after application 
of the scaling factors the Monte Carlo simulations provide a good description of the 
measured data. This allows to determine the set of discriminant variables values which 
provide the best signal to background ratio. 

As a next step we optimize the event selection in order to reduce the background as 
strongly as possible while keeping high signal efficiency. To this end the following cuts 
were varied: 

• the of the kinematical fit 

• the topological AE/aE cut 

• signal box definition in the (xL? xln) plane 

• Rfnin • 
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Figure 6.10: The inclusive distributions of Rmin before (a ) and after (b ) correction. 
Plots c and d present /3cr spectrum for events with charged particles before and after 
the correction, respectively. The red points represent the data, while blue histograms 
denote simulation results. 
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Each set of the cut values was applied to the simulated background and signal sam- 
ples excluding events with charged particles coming from the vicinity of the inter- 
action region and for several sets of tight cuts on the reconstructed Kl energy and 
velocity. This procedure allowed to determine the number of selected background 
events B and the signal efficiency e^T^ as a function of the five variables listed above: 
B = ^E/aE, xL' xL, Rmin) and eg^ = e3^(x^ ^E/aE, xL' xL. Rmin)- Since 

the expected number of events at the end of the analysis chain is small we define the 
following function: 

tcut = , (6.11) 

where Nup denotes the mean upper limit on the expected number of signal events 
calculated at 90% confidence level assuming well - known number of background 
events [53]. The best choice of cut values is defined as the one which minimizes 
the fcut value. As the result of the optimization we have obtained the following values 
of discriminant variables: 

< 57.2 
AE/aE > 1.88 
4.0 < xl < 84.9 (6.12) 
xL < 5.2 

R 

rmn 

> 65 . 

The signal efficiency corresponding to this set of cuts amounts to: 

ea. = 0.233 ± 0.012,^^. 



6.4 Counting of the Ks — > Stt^ events 

After validation of the Monte Carlo simulations and determination of the optimal set of 
cut values defined in Eq. 6. 12 we preform the discriminant analysis of the experimental 
six gamma sample preselected using the requirements for the Kl mentioned before: 

Ecr > 150 MeV 
0.200 < (3^r < 0.225 , (6.13) 

The same analysis was also applied to the simulated background events with six re- 
constructed gamma quanta fulfilling the tight requirements for listed above. This 



''The detailed description of the meaning of the mean upper limit and the statistical methods for its 
estimation will be described in chapter 8. 
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Figure 6.11: Experimental (xL, xL) distribution after the and AE/aE cuts (left). 
The right panel presents the scatter plot of Rmin versus the minimum energy of the 
cluster Emin for events in the Signal Box at the last stage of the analysis. The lines are 
described in the text. 





all 6 - 7 events 


TRV 




AE/a-E 


SBOX 




DATA 


76689 ± 278 


48963 ± 222 


16501 ± 129 


1400 ± 38 


13 ±4 


0± 1 


MC 


76721 ±446 


48984 ± 283 


16230 ± 136 


1210 ±21 


17 ±3 


± 0.06 



Table 6.4: The number of events surviving each subsequent cut. Results for data are 
given in the first row (DATA) and the second row shows the results for Monte Carlo 
simulations (MC). TRV denotes the rejection of events with charged particles com- 
bined with the tight cuts on the energy and velocity of K^. 



provided the estimation of the expected background at the end of the analysis chain. 
The experimental scatter plot of (xin versus xD ^^er the and AE/aE cuts is pre- 
sented in the left panel of Fig. 6.1 1. The solid lines show the signal region defined in 
Eq. 6.12. As it can be seen in the right panel of Fig. 6.11 all events selected by the 
Signal Box are characterized by Rmin less than 65 cm denoted by the dashed line, thus 
at the end of the analysis chain we have found = of the Ks — Stt^ candidates in 
data. The expected background contribution amounts to B^xp = 0. 
For the final cross-check of the simulations credibility at each stage of the selection 
the number of surviving events of both data and Monte Carlo samples were counted. 
These numbers are reported in Tab 6.4. The statistical uncertainties for the results of 
simulations were estimated taking into account the scaling factors and using the for- 
mula presented in Eq. 6.10. It can be seen that the simulated background is consistent 
with data after each cut. Moreover, the agreement has been found also both at the be- 
ginning and at the end of the analysis in all control boxes in the (xItd xD distribution 
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SBOX 


DOWN 


UP 


CDOWN 


CUP 


CSBOX 


DATA 


220 ± 15 


137 ± 12 


5 ±3 


6931 ± 83 


15179 ± 123 


26491 ± 163 


MC 


232 ± 11 


100 ±7 


4±3 


6797 ± 76 


14906 ± 116 


26962 ± 169 



Table 6.5: The population of control boxes in the (xln^ xL) plane right after rejection 
of events with charged particles imposing the tight cuts on the reconstructed energy 
and velocity. 





SBOX 


DOWN 


UP 


CDOWN 


CUP 


CSBOX 


DATA 


13 ±4 


0± 1 


Oil 


Oil 


Oil 


1387 i 37 


MC 


17 ±3 


± 0.06 


± 0.06 


± 0.06 


i 0.06 


1194 i 21 



Table 6.6: Population of control boxes in the (xln^ xD pl^ne defined in Subsec. 6.2.1 
before the cut on Rmin- 

(see Tab. 6.5 and Tab. 6.6). 

Since with the optimal cuts no events were observed only the upper limit on the 
Ks — 37r° branching ratio can be determined. The procedure used to estimate the up- 
per limit taking into account the statistical and systematical uncertainties is described 
in chapter 8. 
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Chapter 7 



Systematic error estimation 

In this chapter we present evaluation of the systematic uncertainties for the measure- 
ment of the Ks — )■ 37r° branching ratio. They are related mainly to the determination 
of the selection efficiencies for the signal and normalization sample €2-^ and es^r, and 
estimation of the background and cuts used in the discriminant analysis. Moreover we 
discuss small corrections due to the differences in the efficiencies of the Ks tagging 
and preselection with so called FILFO filter for the Ks — 37r° and Ks 27r° de- 
cays. FILFO {FILtro FOndo: background filter) is an off - line procedure identifying 
background events at a very early stage of the data reconstruction using only informa- 
tion from the calorimeter. Events rejected by FILFO do not enter the track fitting and 
pattern recognition algorithms which saves CPU-time during events reconstruction or 
reprocessing [47]. The efficiency of the trigger and cosmic veto for both channels has 
been neglected due to the fact, that they were found to be very close to 100% in previ- 
ous KLOE analysis^ [16]. We conclude giving a summary of the estimated systematic 
error affecting our result. 

7.1 Systematics related to acceptance for the Ks — > 27r° 
channel 

For the normalization sample we have considered the following list of systematic ef- 
fects related to the determination of the e2-K efficiency: 

• Difference in splitting and accidental probabilities between data and simu- 
lations 

As it was mentioned in chapter 5, the probability to find one or more accidental 

' Since we have used the tighter „Kl - crash" requirements the energy released in the calorimeter 
by Kl meson was much larger thus the efficiency of trigger is even higher than in the prior analysis. 
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clusters in the acceptance is slightly different for data and simulations. To esti- 
mate the systematic uncertainty originating from this discrepancy we consider 
the fractions of reconstructed number of photons for the simulated Ks — t- 2tx^ 
events -F^™^ = gathered in Tab. 5.2 as the true ones. The presence of addi- 
tional accidental clusters changes the true fraction F^^^ (neglecting the second 
order effects) to: 

F'^^ ^ F^;' ■ (1 - Pai) + ^3*7'^ ■ Pai, (7.1) 

where Pai denotes the probability to find one accidental cluster in the event 
(see Tab. 5.1). This corresponds to the change AF^^ = (Fg*;"^ - F^™^) ■ Pal 
Therefore the systematic uncertainty originating from the difference of the prob- 
abilities APai for data and Monte Carlo amounts to: 

Ae2 Je2. = APai ■ (F*;"^ - Fi;"^)/^!^^ = 7 ■ lO'l 

Similarly we have calculated the systematic uncertainty due to different proba- 
bilities of shower fragmentation for data and simulations. In this case the varia- 
tion of the true F4™^ amounts to: 

AF4^ ^ (3 ■ Fg^^ - 4 ■ Fi;"'^) ■ Psi , (7.2) 

and the corresponding systematic uncertainty: 

A62je2. = APsi ■ (3 ■ F*;- - 4 ■ Fi:;^)/Fi:;^ = 3 ■ lo-^ 

where AP51 is the corresponding difference in the probabilities to generate one 
splitted cluster for data and Monte Carlo. 

• Correction of the cluster reconstruction efficiency 

The systematic error due to this correction was estimated conservatively as the 
difference between the true f^^ fractions evaluated with and without correction 
and amounts to: Ae27r/e27r = 5.2 ■ 10^^. 

• Acceptance related effects 

The number of events counted as a normalization sample was determined in 
chapter 5 taking into account only events with four reconstructed gamma quanta 
Nnorm = N2-k/^2-k- Assumiug, that events with 3-6 reconstructed photons orig- 
inate from the Ks 27r° decay^ Nnorm should be consistent with the result 
determined based on the number of events with 3-6 reconstructed photons: 

^ This assumption was checked looking to the true event decay chains for the simulations. 
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k 





1 


2 


3 


4 


5 


6 


y6 




141 


337 


2916 


16096 


40048 


649 


39 


56832 




29 


43 


352 


352 


148 


21 


1 


522 




83 


89 


89.2 


97.86 


99.63 


96.9 


98 


99.09 



Table 7.1: The number of simulated Ks ~^ 27r° events accepted (N^) and rejected 
(A^^) by the FILFO filter as a function of the number of reconstructed gamma quanta 
k. The numbers were obtained without any cuts on the Kl energy and velocity. 



Source 




Accidental 


0.07 


Splitting 


0.03 


Accept, rel. 


1.50 


Clu. eff. corr 


0.52 


FILFO 


0.46 


TOTAL 


1.65 



Table 7.2: Summary table of the systematic uncertainties on the selection efficiency 
for the Ks 27r° normalization sample. 

^'norra = ^3-67/e3-67^ • Thus, the difference between these two numbers consti- 
tute the measure of the systematic error amounting to Ae27r/e27r = 1-5 ■ 10~^. 

• FILFO preselection efficiency 

To reject the DA$NE background before the track reconstruction a fast fil- 
ter FILFO [47] based only on information from calorimeter is applied. Fur- 
ther reconstruction is done only for events which pass this filter. In the Monte 
Carlo we keep however all the simulated events which allows us to estimate 
the efficiency. To this end we have considered the sample of Ks — 27r° events 
simulated without „Ki - crash" requirements retained and rejected by FILFO 
(see Tab. 7.1). The efficiency is defined as a ratio of the number of events 
accepted by the filter with photon multiplicities 3 - 6 to the total number of 
events with gamma quanta in the same range of multiplicity"^ (see Tab. 7.1), 
and amounts to: = 0.9909 ± 0.0004. This value will be used in the fi- 
nal evaluation of the upper limit presented in chapter 8. As a systematic er- 

^ The £3-67 was determined using the same simulated Ks — > Stt*^ sample which was used for the 
e27r determination (see Tab. 5.2). 

"'We do not consider events with the multiplicities less than 3 because they originate mainly from the 
wrong To time assignment to the event. 
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ror contribution related to the preselection with FILFO we take conservatively 



Summary of the different contributions to the systematic uncertainty on is pre- 
sented in Tab. 7.2, where the total error was evaluated adding all the contributions in 
quadrature. 

Finally, the estimated selection efficiency for the normalization sample amounts to: 



For the search of the Ks — > Svr^ signal the main sources of systematic uncertainties 
originate from the estimation of background and selection efficiency. As in the case of 
the Ks Stt*^ channel we have considered also the systematic effects related to the 
acceptance. 

For systematic study of the background we have repeated the analysis changing the 
parameters values used in xI-k xL calculation (see Tab. 6.1) as well as varying all 
the corrections applied to the Monte Carlo simulations, namely: 

• using the same resolutions for data and simulations in the xln definition 

• using different resolutions for data and simulations in the xlw definition obtained 
with control sample consisting of events with charged particles (mainly Ks — >■ 

tt+7t~;Kl 37r°) 

• removing correction on Rmin 

• repeating the analysis with different energy scale corrections 

• removing correction on the rejection of events with charged particles 

• varying as in the AE/as definition (aE ± S{aE), where 5{aE) denotes the 
standard deviation of as) 

• varying the cuts on the reconstructed K^ energy Ecr and velocity (3cr arbitrarly 



The full analysis was repeated in total twenty times performing each time one of 
the systematical checks listed above. For all of the checks we have not observed any 
changes in the number of background events at the end of the analysis chain. 
As in the case of the Ks 27r° decay we have considered the following systematic 
effects for the selection efficiency related to the acceptance cuts: 



Ae2./e2. = (1 - el)/2e2n = 4.6 ■ lO'l 



= (0.660 ± 0m2stat ± 0.010,^,) . 



(7.3) 




by ± 5%. 
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k 





1 


2 


3 


4 


5 


6 


7 


8 




2 


3 


1 


12 


133 


435 


597 


11 


1 




0.17 


0.25 


0.08 


1.0 


11.1 


36.4 


50.0 


0.92 


0.09 



Table 7.3: The number of events NJl^'^ reconstructed with a multiplicity of clusters k 
for Ntot = 1 195 of KsKl G-fKi events simulated with 0.200 < < 0.225 and 
Ecr > 150 MeV. F*™^ denotes the true fraction defined as F*™^ = iV[f7 ^^^g Nj:^". 



• Splitting and accidental probabilities for data and simulations 

Based on the same simulated sample of the Ks — )■ 37r° events which has been 
used for estimation of the selection efficiency es^r we have determined the true 
fractions -F^™^ of the reconstructed number of photons for signal (see Tab. 7.3). 
The changes of the true Fg™*^ fraction introduced by the presence of accidental 
clusters amounts to approximately: AFq^ = (Fg™^ — Fg™*^) -Pai- Systematic un- 
certainty originating form the difference of the probabilities for data and Monte 
Carlo simulations amounts to: 

Similarly the systematic uncertainty corresponding to the difference in the split- 
ting probabilities has been found to be: 

Aeajeg. = AP^i • (5 • F*;- - 6 ■ Fg*;-)/Fg*;- = 3 ■ 10"^ 

• Correction of the cluster reconstruction efficiency 

The systematic error due to this correction was estimated as in the case of the 
Ks — )■ 27r° channel and amounts to Ae^j^/es-K = 2 ■ 10"'^. 

• Energy scale correction 

As it was described in chapter 5 we have modified the energy scale of the re- 
constructed gamma quanta in the Monte Carlo simulations to provide a better 
agreement with data. Systematic error connected with this correction has been 
estimated by changing the value of the energy shift from 2.2% to 2.6%, and 
amounts to Ae^^/e^T, = 1 ■ 10"^. 

• FILFO preselection efficiency for the Ks — > 37r° 

The efficiency of FILFO filter was estimated analogously to the derivation pre- 
sented in the previous section. The distribution of the Ks Sn^ events sim- 
ulated without any requirements for K^ energy and velocity surviving and re- 
jected by FILFO is presented in Tab. 7.4. The estimated efficiency of the fil- 
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k 





1 


2 


3 


4 


5 


6 


7 


y6 


NX 


11 


12 


1 


23 


208 


658 


903 


15 


1792 
















6 


4 


1 


1 


11 


NX/{NX + Nj^)[%] 


100 


100 


100 


100 


97.2 


99.4 


99.9 


93.8 


99.4 



Table 7.4: The number of simulated Ks 37r° events accepted (A^^) and rejected 
(A'^^) by the FILFO filter as a function of the number of reconstructed gamma quanta 
k. The numbers were obtained without any cuts on the energy and velocity. 

ter for signal is equal to efj^ = 0.994 ± 0.002. This value will be used in the 
final evaluation of the upper limit presented in chapter 8. As a systematic er- 
ror contribution related to the preselection with FILFO we take conservatively 

AeaJes. = (1 - efj/2e3. = 3 ■ IQ-^. 

The last group of systematic uncertainties is connected with the cut sequence used in 
the discriminant analysis: 

• Energy resolution 

The systematical uncertainty due to energy resolution was determined by es- 
timation of the selection efficiency 63^ with different aE values in the AE/aE 
definition, which was varied as it was described in the case of systematics related 
to the background estimation. In this case it amounts to Aes^/ea^ = 1.1 ■ 10~^. 

* -^min cut 

The systematic effects related to the Rmin cut was studied comparing the selec- 
tion efficiencies evaluated with and without the Rmin correction. Difference of 
these two values gives the systematic uncertainty equal to Ae^^^/e^^ = 9 ■ 10""^. 

Systematic effects due to the xju '^^t were investigated using the Ks Stt*^ 
events with four reconstructed gamma quanta. Since the photon multiplicity and 
energy spectrum of the Ks — J- 27r° events differs from the ones for signal we 
expect differences in the shape of x^j^^/ndop for the two samples. To estimate 
the systematic uncertainty related to this difference we have compared the sim- 
ulated x}it/''^dof distributions for the Ks — )■ 27r° and Ks — )■ Sir^ events. For 
both distributions the cumulative curves /g/'" and /g^*" were determined (see 
Fig. 7.1a). The ratio fin'^/fin'^ which is shown in Fig. 7.1b constitutes the esti- 
mation of the contribution to the systematic error. For the cut value used in the 

^ jidof denotes the number of degrees of freedom which amounts to 1 1 for the Ks Stt" events 
and 9 for the Ks —?■ 27r" channel. 



7.2 Systematics related to the selection efficiency and background for the 

Ks 37r° channel 61 




1 2 3 4 5 



9 10 0123456789 10 



%'/ndf 




1 2 3 4 5 



9 10 12 3 



%'/ndf 



x'/ndf 



d) 



5 6 7 8 9 10 

x'/ndf 



Figure 7.1: a) The cumulative curves of the 'x^f^^/ndof distributions for the simulated 
Ks 37r° (black squares) and Ks 27r° (blue circles) samples and b) their ratio; 
c) Comparison of x^^^^/ndof cumulative curves obtained with data (red squares) and 
simulations (blue circles) for the four gamma events and d) their ratio. 



analysis {x^j^^/ndof = 5.2) it corresponds to 1.22 ■ 10"^. Also for the measured 
Ks —7- 27r° events we have constructed a cumulative curve Z^''*" presented in 
Fig. 7.1c. The ratio /i^'**"//^^'^ gives us the second part of the systematic error 
equal to 8 ■ 10~^ (see Fig. 7.1 d ). Adding the two contributions in quadrature we 
obtain the total systematic uncertainty on the selection efficiency related to the 
xjit cut equal to Aes^/es^ = 1.46 ■ 10"^. 



Different contributions to the systematic error on e^,, are summarized in Tab. 7.5, where 
the total systematic uncertainty was evaluated adding all of them in quadrature. 
The final value of the selection efficiency for the Ks Sn^ decay amounts to: 



= 0.233 ± O.OUstat ± 0.006 



(7.4) 
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Source 




Accidental 


0.03 


Splittmg 


0.02 


Energy scale 


l.UU 


\_.iu« cii. curr 


u.zw 


FILFO 


0.30 




1.46 


Energy resolution 


1.10 




0.90 


TOTAL 


2.30 



Table 7.5: Summary table of the systematic uncertainties on the selection efficiency 
for the Ks — )■ 37r° decay. 

7.3 Correction for the different Ks tagging efficiencies 

for the Ks — > 37r° and Ks 2n^ decays 

The difference in the kinematics and in the photon multiplicity between the Ks — )• 37r° 
and Ks — )■ 27r° decays creates a small difference in the Ks tagging efficiencies for 
these channels. This may be result, for example of accidental or splitting clusters which 
can modify spectrum of the reconstructed K^ energy. To take into account this small 
effect we have determined the Ks tagging efficiencies independently for each channel 
using the appropriate simulated sample of events. The estimated efficiencies amounts 
to 

e^; = (23.65 ± 0.12) ■ IQ-^ for the Ks 2n^, 

and el^ = (23.90 ± 0.90) ■ lO^^ for the Ks 3n^. 

The ratio = ^Ir /^Ir = 1.01 ± 0.04 constitutes a correction which will be used in 
the final evaluation of the upper limit on the Ks — )■ 37r° branching ratio. 
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Upper limit on the Kg ^ 
branching ratio and 1 77000 



8.1 Upper limit on the measured number of Ks 37r" 
decays 

The search for the Kg — )• Stt" decay presented in this work failed to detect a signal of 
sufficient statistical significance. Therefore, one can only determine an upper limit on 
the branching ratio for this decay at a chosen confidence level. A limit on a physical 
quantity at a given confidence level is usually set by comparing a number of detected 
events with the expected number of background events in the signal region. The ex- 
pected background depends strongly on the systematic uncertainties existing in the 
measurement, therefore they should be taken into consideration in the limit or confi- 
dence belt calculation [54] . 

In the framework of frequentist statistics confidence limits are set using a Neyman 
construction [55]. This method suffers however from so-called undercoverage when 
the observable is close to the physics boundary (the actual coverage is less than the 
requested confidence level). Moreover, constructed confidence intervals may be un- 
physical or empty, in particular in the case when no events have been observed this 
method gives no answer for the confidence interval [54]. 

Feldman and Cousins [56] proposed a new method of confidence interval estimation 
based on likelihood ratios which automatically provides a central confidence interval 
or an upper confidence limit decided by the observed data itself (so called unified ap- 
proach) [54]. However, if the observable is a Poisson distributed variable, there is a 
background dependence of the upper limit in the case of fewer events observed than 
expected background. This gives a smaller upper limit for measurements with higher 
background, which is clearly undesirable. A solution to this problem was proposed by 
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Roe and Woodroofe taking advantage of a fact that given an observation n, the back- 
ground b cannot be larger than n [57]. Therefore, the usual Poisson pdf (probability 
density function) used to construct the confidence interval should be replaced by a 
conditional pdf. This approach solves the background dependence of the upper limit, 
however, does not satisfy all the requirements of proper coverage and has problems 
when applied to the case of a Gaussian distribution with boundaries [54]. 
All approaches described above do not take into account the systematic uncertainties 
of estimated background and signal efficiency. Several methods have been developed 
to incorporate the systematic errors to the calculation of upper limits. An entirely fre- 
quentist approach has been proposed for the uncertainty in the background rate predic- 
tion [58]. It is based on a two-dimensional confidence belt construction and likelihood 
ratio hypothesis testing. This method treats the uncertainty in the background as a sta- 
tistical uncertainty rather than as a systematic one [54]. 

Several methods combine classical and Bayesian elements (so called semi-Bayesian 
approaches) for example incorporating systematic uncertainty by performing average 
over the probability of the detection efficiency [59]. This method is however of limited 
accuracy in the limit of high relative systematic uncertainties [54]. Conrad et al. ex- 
tended the method of confidence belt construction proposed by Feldman and Cousins. 
This approach includes systematic uncertainties of both the signal efficiency and back- 
ground prediction by integration over a pdf parameterizing our knowledge about the 
sources of the uncertainties [60]. 

In the framework of Bayes statistics the systematic error is included by modification 
of the usual conditional pdf with additional probability density functions for the back- 
ground expectation and signal efficiency [54, 61]. 

Summarizing, there are many approaches for upper limit estimation and still there is a 
lot of effort made towards improving these methods and understanding their practical 
implications. All the approaches mentioned before give in principle different confi- 
dence intervals and one has to choose the method depending on the relevance of sys- 
tematic uncertainties and personal opinion about the Bayesian and classical statistics. 
The analysis described in this thesis resulted in n = observed events with the number 
of expected background events 6 = 0. In this case all the methods reveal that the upper 
limit is almost insensitive to the systematic errors of the background estimation and 
signal efficiency [54, 58, 59, 61]. Moreover, in the case of n = 6 = the upper limit 
on the number of signal events at the 90% confidence level amounts to 2.21 [58] or 
2.33 [54, 59-61] for the systematic uncertainties on background and signal efficiency 
less than 10%. As it was shown in chapter 7 in our search the systematical errors do 
not exceed 5%, therefore we assume the upper limit on the number of the Ks Stt*^ 
events at the 90% confidence level amounting to A''3^(0.9) = 2.33. This number will be 
used in the next section for the calculation of the upper limit on Ks 37r° branching 
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8.2 Determination of the upper limits on the Ks — > Stt 
branching ratio and 1 77000 1 



From the limit on the number of expected signal events one can calculate a limit on 
the branching ratio. Recalling Eq. 5.2 and taking into account the corrections for the 
difference in tagging and FILFO preselection efficiencies for signal and the Ks 2it^ 
normalization sample described in chapter 7 we obtain the following expression for the 
branching ratio: 



BR{Ks ^ 3vr°) = ^ ■ ^ ■ ■ BR{Ks ^ 2vr°) , (8.1) 



where R^r denotes the ratio of tagging efficiencies for signal and the Ks — )■ 27r° 
normalization sample amounting to Rcr = 1-01 ± 0.04. e^^ = 0.994 ± 0.002 and 
ef^ = 0.9869 ± 0.0006 are the FILFO preselection efficiencies determined in chap- 
ter 7. The branching ratio for the Ks 27r° channel is equal to: BR(Ks 27r°) = 
0.3069 ± 0.0005 [19]. Taking into account value of the upper limit on the number of 
expected signal: A^3^(0.9) = 2.33, the selection efficiency for the Ks Sn^ chan- 
nel: 63^ = 0.233 ± 0.012stat ± O.OOSsys and the total number of Ks — )■ 27r° events: 
N2n/e2n = (1.142 ± 0.005) ■ 10*^, we have obtained the upper limit on the Ks 37r° 
branching ratio at the 90% confidence level: 



This value is almost five times lower than the latest result BR{Ks — 37r°) < 1.2 ■ 10~^ 
published by KLOE [16]. The upper limit on the Ks — )■ 37r° branching ratio can be 
translated into a limit on the |?7ooo| at the 90% confidence level [37]: 



which corresponds to an improvement of the |?7ooo| uncertainty by a factor of two with 
respect to the latest direct search [16]. 



ratio. 



'cr 
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Chapter 9 

Summary and outlook 



Kaons have been a remarkably important particles in the development of the Standard 
Model. The kaon system, being a relatively simple one, played the key role in the dis- 
covery of such phenomena as parity and CV violation, the GIM mechanism and the 
existence of charm, and was central in the investigation of lepton flavour and CVT 
symmetries [62]. 

CV violation is deeply related to such fundamental issues as the microscopic time re- 
versibility of physical laws and the origin of the baryonic asymmetry of the universe. 
It is also the only known phenomenon which allows an absolute distinction between 
particles and antiparticles [63]. Since the first discovery in the neutral kaon system the 
CV symmetry breaking has been a very active field of research. Although at present 
the main experimental effort is focused on the neutral B and D meson system studies, 
there are still several interesting open issues in the kaon physics demanding investi- 
gation. One of them is the Ks — > 37r° decay which still remains undiscovered. The 
best published upper limit on its branching ratio BR{Ks — > 37r°) < 1.2 ■ 10"^ is 
still two orders of magnitude larger than the predictions based on the Standard Model. 
Therefore, the complete understanding of the CV violation in the neutral kaon system 
demands a new high - precision experiments which will contribute also to the CVT 
conservation tests. 

This work presents the search of the Ks — 37r° decay based on the data sample gath- 
ered in 2004 - 2005 with the KLOE detector operating at the $ - factory DA$NE in 
the Italian National Center for Nuclear Physics in Frascati. DA<I>NE is an e+ and e~ 
collider optimized to work at the center of mass energy ^/s = 1019.45 MeV. In the two 
storage rings of DA$NE 120 bunches of electrons and positrons are stored. Each bunch 
collides with its counterpart once per turn, minimizing the mutual perturbations of col- 
liding beams. The e'^e~ collisions result in the meson creation which is almost at rest 
and decay predominantly to kaon pairs. The decay products are registered using the 
KLOE detection setup surrounding the e+e" interaction point. KLOE consists of large 
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cylindrical drift chamber surrounded by the electromagnetic calorimeter. The detec- 
tors are immersed in the axial magnetic field generated by superconducting solenoid. 
The Ks mesons were identified via registration of these mesons which crossed the 
drift chamber without decaying and then interacted with the KLOE electromagnetic 
calorimeter. The Ks four - momentum vector was then determined using the regis- 
tered position of the Kl meson and the known momentum of the meson. The search 
for the Ks — )■ 37r° — 67 decay was then carried out by the selection of events with six 
gamma quanta which momenta were reconstructed using time and energy measured by 
the electromagnetic calorimeter. Background for the searched decay originated mainly 
from the Ks — 27r° events with two spurious clusters from fragmentation of the elec- 
tromagnetic showers or accidental coincidences with signals generated due to particle 
loss of DA$NE beams, or from false K^ identification. To increase the signal over 
background ratio after identification of the Ks meson and requiring six reconstructed 
photons a discriminant analysis was performed. It started from rejection of events with 
charged particles coming from the vicinity of the interaction region which suppress the 
(f) — )■ KsKl — )■ {Ks — )■ ii^Ti:~ ,Kl — Svr^) background events. The further analysis 
was based on kinematical fit, testing of the signal and background hypotheses and ex- 
ploiting of the differences in kinematics of the Ks decays into 2tx^ and Svr^ states. 
As a result of the conducted analysis no events corresponding to the Ks — t- Stt^ de- 
cay have been identified. Hence, we have obtained the upper limit on the Ks Stt*^ 
branching ratio at the 90% confidence level: 



This value is almost five times lower than the latest published result. This upper limit 
can be translated into a limit on the |?7ooo| the 90% confidence level: 



which corresponds to an improvement of the |?7ooo| uncertainty by a factor of two with 
respect to the latest direct measurement [16]. 

However, the upper limit on the Ks — )■ Stt" branching ratio determined in this work 
is still about one order of magnitude larger than the prediction based on the Standard 
Model. Thus, the picture of CV symmetry violation in the neutral kaon system re- 
mains incomplete. Therefore, among several other experiments aiming in the precise 
measurements of rare and ultra - rare kaon decays [62], we are continuing the research 
of the Ks 37r° decay by means of the KLOE-2 detector. The KLOE-2 collaboration 
is continuing the physics program of its predecessor. In the last years, a new machine 
scheme based on the Crab - waist optics and a large Piwinsky angle [64] has been pro- 
posed and tested to improve the DA$NE luminosity. The test has been successful and 
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presently DA$NE can reach a peak luminosity of a factor of three larger than previ- 
ously obtained. The next data taking campaign during 2013-2015 will be conducted 
with a goal to collect total integrated luminosity amounting to about 20 fb~^, which 
corresponds to one order of magnitude higher statistics with respect to what was used 
in this search. For the forthcoming run the KLOE performance has been improved by 
adding new subdetector systems: the tagger system for studies of the meson produc- 
tion in the 77 reactions, the Inner Tracker based on the Cylindrical GEM technology 
and two calorimeters in the final focusing region [65, 66]. These new calorimeters will 
increase the acceptance of the detector, while the new inner detector for the determi- 
nation of the Ks vertex will significantly reduce the contribution of the background 
processes involving charged particles. Increasing the statistics and acceptance of the 
detector while significantly reducing the background gives the realistic chances to ob- 
serve the Ks 37r° decay for the first time in the near future. 
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